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Abstract 


In this thesis, the linear transport equation with the entry condition given on the entry part of 
the boundary is considered and an estimate of the solution in WL is proved. Further, by using 
this result, the existence and the uniqueness of the local solution to a quasi-linear equation 
system arising from the modeling of the atmosphere with phase transition of water are proved. 
Moreover, the possibility of weakening the condition on the velocity field on the entry part of 


the boundary for obtaining an estimate similar to that of the first result is shown. 


Key-words: Transport equation, entry condition, method of characteristics, model of atmo- 


spher. 


ii 


Résumé 


Dans la présente thèse, on considère l’équation de transport linéaire (équation aux dérivées 
partielles du premier ordre de type hyperbolique) avec la condition d’entrée donnée sur la par- 
tie entrante de la frontière et on démontre une estimation de la solution dans W}. De plus, 
en utilisant ce résultat, on démontre l’existence et l’unicité de la solution locale d’un système 
d'équations quasi-linéaire issu d’un modèle décrivant le mouvement de l’air avec la transition 
de phase de l’eau. En outre, on affaiblit la condition sur le champ de vitesse sur la partie en- 


trante de la frontière, pour l'obtention d’une estimation similaire à celle du premier résultat. 


Mots-Clés : Équation de transport, condition d’entrée, méthode des caractéristiques, modèle 


pour l'atmosphère. 
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Introduction 


Notes historiques sur l’équation de transport 


dl N 1890, en se basant sur la théorie développée par Lagrange et complétée par 
Cauchy et par Jacobi, Edouard Goursat a élaboré dans son ouvrage [GOU 21], 
une méthode à l’aide de laquelle on ramène l'intégration d’une équation aux dérivées 
partielles du premier ordre à l'intégration d’un système d'équations différetielles or- 


dinaires à une seule variable indépendante. En effet, on considère en particulier une 


équation à deux variables indépendantes 


PE. () 
"ôx 0y 


où les coefficients A1, 4, B ne dépendent que des variables x et y et de la fonction 
inconnue z. L'équation (1) traduit le fait que le plan tangent à une surface intégrale 
Z(x, y) en un point (x, y, z) de cette surface contient la droite (D), passant par ce même 


point, et qui a pour équations 


Selon Goursat, une courbe caractéristique est une courbe de l’espace qui, en chacun 


de leurs points, est tangente à la droite (D) correspondante. De cette définition on peut 
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tirer la propriété :“ Toute surface intégrale est un lieu des courbes caractéristiques et in- 
versement, toute surface engendrée par des courbes caractéristiques est une surface in- 
tégrale”. Par conséquent, il a constaté alors que sur une caractéristique, les paramètres 
directeurs de la tangente dx, dy,dz, sont proportionnels aux coefficients A1, 42,B et 


ces courbes sont donc définies par le système d'équations différentielles ordinaires 


dx _ dy _dz 


A) A  B 

Plus tard, en 1937, Courant et Hilbert ont consacré un chapitre tout entier dans 
leur fameux livre [COU 89] à la théorie des équations aux dérivées partielles du premier 
ordre où ils donnent une description, de la méthode des caractéristiques, simillaire à 
celle de Goursat. Dans la période entre 1939 et 1969, les problèmes de transport ont été 
largement étudiés et cela bien sûr est motivé par les besoins de la technologie. C’est 
pourquoi l’équation de transport linéaire, à cette époque, paraît dans des domaines 
très variés, elle provient principalement de la physique en particulier dans la neutro- 
nique ([DAV 57], [OR 58], [REE 65], [ZWE 67]), le transfert radiatif ([CHA 50], [STE 671), 
la physique des plasmas ([MON 64], [OBE 67]) et la mécanique des fluides ([SER 591], 
[HAR 67]). 

Dès lors, un grand nombre de mathématiciens s’intéressait à l’étude de l’équa- 
tion de transport dans ses différentes formes et ses multiples applications, citons par 
exemple [BAR 70], [BLO 70], [Drp 89], [AMB 04], [BOY 05], [CRI 14], [ASC 14]. 

Dans son travail [BAR 70], Bardos a examiné l'équation de transport avec la condi- 
tion d’entrée sous des conditions suffisamment générales mais avec le champ de vi- 
tesse de transport v appartenant à la classe €! ; il utilise en substance la méthode des 
caractéristiques, mais avec des "méthodes modernes"-analyse fonctionnelle et théorie 
des semi-groupes-. 

Dans le cas où v est moins régulier (appartient à un espace de Sobolev par exemple), 
Diperna et Lions [DIP 89] ont introduit la théorie des solutions renormalisées à l’aide 


de laquelle ils ont démontré l’existence et l’unicité des solutions faibles dans L® d’une 
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équation de transport de la forme suivante 
Ô;u—v-Vu+cu=0 (2) 


avec une donnée initiale dans L®. En effet, ils ont établi que toutes les solutions faibles 
de (2) sont des solutions renormalisées. Plus précisément, ils ont énoncé que pour 


toute fonction régulière B, la fonction B(u) est une solution de l’équation de transport 
04B{u) — v-VB(u) + cuB'(u) = 0. 


Cette théorie a été étendue au cas d’un champ de vitesse appartient à l’espace des fonc- 
tions à variation totale bornée (BV) par Ambrosio [AMB 04], où il a traité la question 
de l’unicité des solutions faibles bornées du transport par un champ BV. On cite éga- 
lement le travail [CRI 14], dans le quel les auteurs ont étudié un problème aux limites 
pour une équation de continuité dans le cas d’un champ de vitesse borné, à divergence 


bornée et est dans L!_ (0, T; BV). 


loc 


Motivation physique et mathématique 


La loi de la conservation, une des lois fondamentales de la mécanique des fluides, 


est exprimé par l'équation dite équation de continuité 
06 +V-(pv) =0, (3) 


où p est la densité du fluide en considération, tandis que v = (V1, v2, va) est la vitesse 
avec laquelle le fluide se déplace. Cette équation peut être considérée comme proto- 
type des équations de transport, en particulier de celles qui se considèrent dans un 
milieu fluide. En effet, si on considère une substance diffusée dans un fluide qui se dé- 


place suivant son mouvement, comme une substance dissoute dans le fluide, la den- 
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sité de cette substance w vérifierait une équation du même type que (3) 


OF + V-(yv) = 0. 


Or, cette substance peut subir des réactions chimiques ou d'autre processus, qui peuvent 
provoquer la variation de la masse de cette substance. Donc souvent on devra consi- 


dérer une équation du type 


dy +V-(yv) = f(Y,x, LE} (4) 


où f(Y, x, t) est le terme de la variation de densité due au processus complémentaire. 
Naturellement, pour chaque substance et pour chaque problème spécifique, on devra 
préciser la forme de la fonction f(w, x, ©). 

Par exemple, si on considère ce qui se passe dans l’atmosphère, on voit que la va- 
peur d’eau et les gouttelettes d’eau joueront le rôle principal pour la formation des 
nuages et la précipitation. La vapeur d’eau se déplace avec l'air, tandis que les goutte- 
lettes d’eau se déplacent sous l'influence du mouvement de l’air ainsi que par la force 
gravitationnelle. La vapeur d’eau peut se transformer en eau liquide par le processus 
de condensation et former des gouttelettes et les gouttelettes d’eau peuvent être trans- 
formées en vapeur d’eau par l’évaporation. En outre, les gouttelettes peuvent subir le 
processus de coagulation entre elles. Ces circonstances nous conduisent à des équa- 
tions du type (4) avec le terme de processus complémentaire f(w, x, f) bien spécifique. 

À titre d'exemple, nous citons la modélisation des phénomènes atmosphériques 
proposée dans [SEL 11]; nous allons aussi appliquer notre résultat “théorique” à ce 
modèle. Dans le modèle proposé dans [SEL 11], les auteurs ont cherché à décrire le 
mouvement de l’atmosphère avec la transition de phase de l’eau. On considère en effet 
un système de trois équations de conservation de la masse (de la forme (4)) pour: l’air 
sec, la vapeur d’eau et l’eau liquide contenue dans les gouttelettes. Par “air sec” on 


entend la partie de l’air consistant dans toutes les composantes sauf H,0; c'est-à-dire, 
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la partie composée par No, O>, Ar, CO, etc. mais sans considérer H,0. Pour l’air sec, la 
loi de la conservation de la masse est exprimée par l'équation de continuité (3), c’est- 
à-dire 


ô6+V-(pv) =0, (5) 


où p est la densité de l’air sec, tandis que v désigne la vitesse de l’air. Pour la densité 
7x de vapeur d’eau, si nous désignons par H,, la quantité (par l'unité de volume et par 
l'unité de temps) de H,0 qui se transforme du gaz au liquide et par H,, la quantité qui 


se transforme du gaz au solide, on aura l'équation 
On +V:(rv) = - Hg] - Hgs. (6) 


Si la température est supérieure à la température de fusion (T > 0°), on a H,, = 0, de 


sorte que l'équation (6) se réduit à 
On +V-(rv) = -H](r, 0). (7) 


D'autre part, en désignant par w la vitesse des gouttelettes d’eau et en considérant la 
condensation ou l’évaporation produite sur les gouttelettes représentée par la fonction 
hgi, on donne l’équation pour la densité o de l’eau liquide contenue dans les goutte- 
lettes de masse m 


00 +V-(OW) +0m(Mmhg10) — hg1O = (8) 
=-0G\(7; m) - m [ Btm,m')o(m)o(m')dm'+ 
0 
+Go(x, 0; M) + = | Bim-m',mo(m')o(m-m'dm', 
0 


où G1(7; m) représente le taux de disparition de gouttelettes de masse m, Go(x, 0; m) 
est le taux de création de gouttelettes de masse m et B(m, m') est la probabilité de 
rencontre entre une gouttelette de masse m et une de masse m’. 


Du point de vue mathématique, si le champ de vitesse v est donné et jouit de la 
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régularité suffisante, les équations du type (4) avec un f(y, x, t) suffisamment simple 
peuvent être résolues par la méthode classique des caractéristiques. C'est-à-dire, en 


définissant les caractéristiques y par l'équation 
d 
—Y(#) = v(E,Y(D)), 
Fe mi At V(E, (0) 


le long de chaque caractéristique y, on transforme l'équation (4) en une équation dif- 


férentielle ordinaire 


d 
Ve up + ft, 7,0. (9) 


La résolution de (9) donnera la résolution de (4). 

Toutefois, dans les applications d'intérêt physique, on rencontre souvent le cas où 
la régularité de v (et des autres termes) est insuffisante pour appliquer directement la 
méthode des caractéristiques. Ainsi, nous sommes obligés à recourir à une construc- 
tion de la solution plus complexe, pour laquelle l’étude de l’équation linéaire et l’es- 
timation de sa solution sont fondamentales. Par exemple, la résolution du système 
d'équations non linéaires (5)-(8) se base sur l'étude des équations linéarisées, en par- 
ticulier, les estimations utiles des solutions des équations linéarisées. Nous rappelons 
aussi que pour l'application aux équations non linéaires, on trouve souvent particuliè- 
rement utiles les estimations du gradient de la solution dans L®. 

L'autre aspect important est la condition d'entrée. En effet, l'hypothèse de l’annu- 
lation de la composante normale de la vitesse v sur la frontière du domaine dans le- 
quel on considère l'équation facilite le calcul et le raisonnement, comme les auteurs 
de [ASC 14] l’ont fait dans leur travail. La condition que la composante normale de la 
vitesse de transport sur la frontière du domaine considéré n’est pas nécessairement 
nulle, c’est-à-dire qu’il y a l'entrée et la sortie de la quantité représentée par la fonction 
inconnue de l’équation, pose des problèmes très délicats sur la condition d’entrée dans 
le voisinage des points où la composante normale de la vitesse de transport s’annule. 


Ces motivations physiques et mathématiques nous ont ramenés à l’étude que nous 
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présentons dans la suite. 


Plan de la thèse 


Le plan de la présente thèse est le suivant : le premier chapitre présente le résul- 
tat de l'estimation dans W} de la solution de l'équation de transport linéaire avec la 
condition d’entrée. Ce résultat (voir [BAZ 18]) généralise le résultat de [ASC 14] au cas 
où la composante normale de la vitesse n’est pas nécessairement nulle. Plus précisé- 
ment, nous allons considérer une équation de transport linéaire avec une condition 
d'entrée sur la partie entrante de la frontière et on établi une estimation de la solu- 
tion dans la classe W1. Comme la condition de régularité du champ de vitesse n’est 
pas suffisante pour utiliser directement la méthode des caractéristiques, nous allons 
construire un problème approché en régularisant les coefficients de l’équation ainsi 
que les données initiale et d'entrée, mais cela se fait d’une manière particulière à cause 
de l'entrée. Ensuite, on établit des estimations des solutions approchées et en passant 
à la limite, on obtient le résultat cherché. 

Dans le chapitre 2, nous présentons une application des résultats obtenus dans le 
chapitre précédent (voir la deuxième partie de [BAZ 181). Plus précisément, en utilisant 
le résultat du cas linéaire, nous allons démontrer l'existence et l’unicité de la solution 
locale d’un système d'équations non linéaires avec l'entrée. C’est un sous-système du 
système d'équations du modèle de l’atmosphère avec la transition de phase de l’eau 
(ce modèle est proposé dans [SEL 11]). 

Dans le chapitre 3, nous proposons d’affaiblir la condition technique posée dans le 
premier chapitre pour l'obtention de l'estimation. En effet, dans le premier résultat, on 
a posé une condition de la régularité du champ de vitesse dans le voisinage de la région 
où sa somposante normale s’annule. Pour affaiblir cette condition, nous proposons la 
régularisation des fonctions données par l’opérateur régularisant variable (qui dépend 
de la position x). Puis, nous démontrons que la mesure de la zone affectée par une ré- 


gularisation particulère tend vers 0, quand le paramètre € de l’opérateur rélgularisant 
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tend vers 0. Ceci nous permet d'obtenir un résultat analogue à celui du premier résul- 
tat sur l'estimation dans W} de la solution de l'équation de transport linéaire avec la 


condition d'entrée. 
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Estimation dans W? de la solution de 
l'équation de transport avec une 


condition d'entrée 


1.1 Introduction 


f E problème de l'existence et de l’unicité ainsi que de la régularité de la solution 
EXD d'une équation de transport linéaire a été étudié par nombre de mathéma- 
ticiens (voir par exemple [BAR 70], [Dip 89], [AMB 04] pour ne citer que quelqu’uns 


parmi les travaux les plus connus). Or, dans le cas où l'équation est envisagée avec 


Chapitre 1. Estimation dans W} de la solution de l’équation de transport … 


une condition d'entrée (donnée à la frontière du domaine), on trouve une probléma- 
tique plus complexe, comme il a été évoqué par plusieurs auteurs [BAR 70], [SEC 98], 
[PLO 12], [CRI 14], [BOY 05]. Dans [PLO 12], les auteurs ont étudié en détail le pro- 
blème pour l'équation de transport stationnaire avec des données d’entrée. En par- 
ticulier, ils considèrent le problème de l’existence et de l’unicité de la solution dans le 
cadre des espaces de Sobolev W;(Q), 1 < p < co, 0 < s< 1; pour éviter l’irrégularité de 
la solution, ils imposent des conditions assez complexes dans le voisinage de la fron- 
tière. Plus récemment, la question de l’affaiblissement de la condition sous laquelle il 
existe une solution a été approfondie de manière fort accentuée. On cite par exemple 
le travail [AMB 04], dans lequel l’auteur considère une équation de transport dont le 
champ de vitesse appartient à l’espace des fonctions à variation totale bornée (BV). 
Plus tard, dans [CRI 14] les auteurs ont établi le caractère bien posé des problèmes aux 
limites pour les équations de continuité où le champ de vitesse est borné, à divergence 


bornée et est dans L},. 


(0, T; BV(Q,;R”)) pour tout ensemble ouvert et borné Q, € Q. 
Or, la question de l’estimation du gradient de la solution n’a - semble-t-il - pas aussi 
attiré l’attention des mathématiciens dans les derniers temps. Mais nous constatons 
que, lorsqu'il s’agit des équations non linéaires, l'estimation du gradient de la solution 
de l’équation linéarisée est souvent la clé de la résolution du problème non linéaire. 
Dans [ASC 14] les auteurs, en étudiant l’équation de transport linéaire, ont établi 
des estimations, y compris l'estimation dans L® du gradient, de la solution. Or, dans 
[ASC 14] on a supposé que la composante normale de la vitesse (du flux representant le 
transport) s’annule sur la frontière du domaine dans lequel l'équation est considérée. 
Dans le présent chapitre on va généraliser le résultat de [ASC 14] sur la norme dans 
L® du gradient de la solution de l'équation de transport linéaire au cas où la compo- 
sante normale de la vitesse n’est pas nécessairement nulle; il s’agit donc du cas où il y 
a l'entrée et la sortie du domaine considéré. Nous allons utiliser les idées de [ASC 141, 
mais la présence de l'entrée nous conduit à une construction particulière d'équations 


approchées. 
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1.2 Position du problème linéaire 


Dans ce chapitre, nous allons considérer une équation de transport linéaire dans 
un sous-ensemble ouvert borné Q de R” (n > 2). Pour la frontière 0Q de Q, nous sup- 
posons qu’elle est de classe €!, Q étant d’un seul côté de 0, et qu’il existe un £0 > 0 


tel que, pour tout x e 00, on ait 
x-Eenñn(x)eQ Veel0,eol, (1.1) 


où ñ(x) désigne la normale extérieure unitaire à la frontière 0Q au point x € 0Q (dans 
la suite nous écrivons simplement n, au lieu de ñ(x)). 


On considère une fonction v définie sur [0,oo[ xQ à valeurs dans R” telle que 


ve L},.(0,00 W, (G;R”)), (1.2) 
V:veL}, (0,00; We (Q)), (1.3) 


où 1<q<o0. 
La fonction v sera considérée comme le champ de vitesse qui définit le flux pour 
une équation de transport. Nous définissons, pour chaque f € [0, col, la partie d’entrée 


F_(®) et la partie de sortie l';(#) de la frontière 0Q par les relations 
F_(1H={xe0dQ]|n,:v(t,x) <0}, (1.4) 


F,(D ={xEe0Q1n,:v(E,x) > 0}. (1.5) 


Dans [0,0co[ x*Q nous allons considérer l’équation de transport linéaire pour une 


fonction inconnue u (à valeurs réelles) 


GER )=gu+ (1.6) 
Gi uv) =gu+f, À 
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où get f sont deux fonctions données à valeurs réelles (v est la fonction introduite 


ci-dessus). L'équation (1.6) sera envisagée avec la condition initiale 


u(0,x) = #p(x) pour xeQ (1.7) 
et la condition d'entrée 
u(t,x)=u(t,x) pour(t,xe [J [{4xr_(|], (1.8) 
0<f<oo 


uo(x) et u1(f, x) étant des fonctions données. On suppose la condition de compatibilité 
Uo(X) = U1(0,x) pour xel-(0). (1.9) 


Nous rappelons que pour le problème (1.6)-(1.8) les points où v-n, s’annule posent 
des questions délicates. Pour cela, nous excluons le cas d’un saut de la valeur de v-:n, 
quand v-n, change de signe et, sur les points (f, x) € [0,co[ x0Q tels que x e 0Q\(T-(H)U 


F}(P)), nous supposons les conditions suivantes : 


CONDITION A. Nous supposons qu'il existe deux constantes £1,€2 > 0, € < 2 (£Eo 


est donné dans (3.1)), telles que v, Dv, Vf, Vg, V(V: uv) soient continues dans X:,, où 
Ze, = {(f, X) € [0,o0[ xQ] dist((#, x), 2) < €2}, (1.10) 


Z={(t,x) € [0,oo[x0Q] —€1 < v(f,x): ny <E1}, (1.11) 
ô 
Dv = Gx V'ji,j=1,..,n 
CONDITION B. Si xo € 0Q\(T_(to) UT + (t)), alors, quelle que soit la caractéristique y 


passant par un point (f, x) € [0,oœ![ xQ, on a 


(Lo; Xo) € Y. 
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Ici la caractéristique y = {(Yo(£),71(#),-+,Yn(0)}rer,, Passant par le point (£, x) € 


[0,oof x Q se définit par le système d'équations différentielles 


0 = : (= v;(yE)), j=1,.. 
ar 70 — d au d = V; Ÿ » J = » »N, 


avec la condition initiale 
Yo =t, YjO=x;, j=1l,...,n; 


OÙ 1nax est l'intervalle maximal d'existence de la solution. 

Pour ce qui concerne le gradient de uw1(f,x) sur l°_(f), outre le gradient Vr (ui 
défini de manière naturelle sur la variété de dimension n — 1, nous défnissons la com- 
posante normale n;:Vu; par 


1 
x Vun|r_ = AU Ur VE U] — ÔU; — UV: V+gU] + f), (1.12) 


où v; =v-(v-n;)n,. La nécessité de la définition (1.12) est motivée par le fait que, si 
on susbstitue v =(v.n,;)n,;+v, dans (1.6),ona 


(v-nx)nx Vu=-v;:Vu-0çu-uV-v+gu+f. 


Dans la suite, pour ne pas alourdir l'écriture, par abus de language, nous utiliserons les 
notations 


IVau Loc = IVr we + Ax: Vuillze(r (0) (1.13) 
Ie lus cp = M llzer-( + IVw lier (y (1.14) 


où ñn,- Vu, est la fonction définie dans (1.12). 
Sur la régularité et éventuellement sur le signe des fonctions f, g, uo, 1, on va les 


préciser dans l'énoncé des résultats. 
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1.3 Résultat principal pour l'équation linéaire 


Le résultat principal de notre travail sur l’équation linéaire (1.6) est le suivant. 


Théorème 1.3.1. On suppose que les conditions (1.2), (1.3) ainsi que les conditions À 


et B sont remplies et que, quel que soit t > 0, les fonctions v, g, f, uo, w1 vérifient les 


relations 

geLT(0,EWL(Q), feLl(0,EWLi(Q), ue WL(O), (1.15) 
v, Vr_çour, f Ou, g—-V:veL®{( (J [{a4xr_(nl]), (1.16) 

O<t<t 
n,-Vaie Le (NOMES net); (17) 

O<t<t 

u,est uniformément continue sur U [{4xT-(#] et (1.18) 

O<t<t 


y admet presque partout la dérivée par rapport à t, 


où1=< q<etn,- Vu, est la composante normale du gradient de u\ dans le sens défini 
par le second membre de (1.12). Alors il existe une fonction u et une seule appartenant 
à la classe 


L°(0, 1; W(Q) n W,(0, 5; L°(Q)), 


qui vérifie l'équation (1.6) presque partout dans [0, t] x Q, la condition (1.7) presque 


partout dans Q et la condition (1.8) presque partout sur U [{4xT-(]. 


O<f<oo 
De plus, la fonction u vérifie l'inégalité 
| U]| 1°0(0,#:w2. (0) <max(A, sup ess B(s)), (1.19) 
O<s<t 
où 
nel 
Aze 7 Vélisorwbont IV loue tlDUliaore an) \ 


q-1 
q 


x [|| Uollw1, (0) +1 Pl Laco,zwz, cop) 
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LE 
B(9) = 60° * UsliczwLoytIVvliaczwt ant lPvliaszreto), 


_ 
x [1241 CS) co + ES 4 face wz cop: 
LT-(S) (SEWI(Q)) 


_ : Et 
dans le cas q = , on remplace = par 1. De plus, si f, uo et u1 sont non-négatives, u 
p a P 


est aussi non-négative. 


On remarque que, à cause du facteur dans le deuxième membre de (1.12), 


Nx 
les conditions (1.16) ne sont pas suffisantes pour garantir l'appartenance du second 
membre de (1.12) à L°(T_(#)) et pour cela nous imposons la condition (1.17). 
Pour démontrer le théorème 1.3.1, nous allons construire les équations approchées 
par la régularisation des coefficients et des données et établir des estimations des so- 
lutions approchées en utilisant les caractéristiques (une méthode dont l’utilisation di- 


recte sur le problème donné est rendue problèmatique par la basse régularité de v); le 


passage à la limite nous donnera le résultat énoncé. 


1.4 Préliminaires 


Pour obtenir des estimations et démontrer l'existence et l’unicité de la solution, 
nous allons utiliser des équations approchées avec les données initiales et d'entrée 


approchées. Pour ce faire, nous écrivons d’abord l'équation (1.6) dans la forme 


dçu+v-Vu+cu=f, c=V:v-g. (1.20) 


Construisons d’abord l’approximation des fonctions v, c, f. 


On commence par prolonger les fonctions v, c, f pour { <0 par 


v(t,:) =0, c(£,:) = C(—f,:), LUS) = 0. 
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On prolonge aussi >, dans [-1,0[ xQ par 


Ze = {(f,X) € [-1,o0[ xQ]dist((#, x), 2) < €2} (1.21) 


avec le même > défini dans (1.11) (et donc 2 € [0,1 x00Q). 


Pour (f,x) ER x R” on pose 


3d(t, x) 
€2 


)), 


h(t, x) = max(0, min(1,2 — 


où d(t, x) = dist{(#, x), ([-1,00[xQ)\3,,) est la distance entre le point (r,x) € R x R" et 


l’ensemble ([-1,0c0[ xQ)\Z,,.On définit 


ht, 2) = (Ce, * h)(,%), (1.22) 


où C4, est un noyau régularisant avec un rayOn €} = 2 (ici la convolution est relative à 
(t, x)). On considère deux fonctions p;(-) et &4(-) (0 < € < min(1, D) définies surRet R” 


respectivement par 


l A l X 
Oe(t) = —-p1(-), ée(x) = oil) (1.23) 
€ € € € 
où 
en 
5 — pour -l1<t<l 
mt=4 Le ras , 
0 ailleurs 
+, 
= pour |x| <1 
1-ly2 q 
é1(X) ES y 
0 ailleurs 


Nous définissons dans R x Q les fonctions 


ve(t,x) = (1 ht, x))u(t, x) + h(t,x)0(1,x), (1.24) 


ce(,x) = (ht, x))c(t, x) + ht, Gt, x), (1.25) 
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ÉGH=QA-hAGO)fE +R fx, (1.26) 


où 
J'éetx- put, y)dy 


De(t,20 =? 


Jau-pay "26 
Q 


J'éetx- pet, dy 


= _Q 
Ce(t, x) = Te *1Qe(0), 
Q 


Jéex-nft dy 


ve Q 

ES e(E), 

FD a — #06 0 
Q 


et +, désigne la convolution par rapport à f. 


Nous définissons aussi dans Q 
uË (x) = (1— h(0,x))uo(x) + A(0, x) (x), (1.27) 


où 
J'éetx- Yuo(y)dy 


ee 
: J&x- dy 
Q 


En ce qui concerne la donnée d’entrée w1, qui est définie sur U [{4xr-(n|], 
0<f<oo 

pour la commodité de notation, nous la prolongeons par 0 sur (f, x) € [0,œ![ x0Q tels 

que x g l'_(f) (en réalité ce prolongement n'intervient pas dans nos raisonnements 


grâce au choix 0 < € < min(I, 2) mais peut faciliter l'écriture); on prolonge u1(f, x) 


aussi sur [—-1,0[x0Q par 0. Cela étant, on définit 


L &e(X— y) u(t,y)dS, 


€ ôQ 
d PE do pds, 
ôQ 


HE (1, x) = 


uË(£, x) = (1 RE, x) (GC, 2 #1 060) (0 + AG, D) GC, 2 +060) (D. (1.28) 


Les fonctions ce, Le, fe, uf, u\ étant bien définies, on considère l'équation pour la 
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fonction inconnue uf 
Ou + ve Vu + cu = fe+uÿpelt) dans [-1,œ[xQ, (1.29) 


avec la condition initiale 


u‘(-1,)=0 sur OQ, (1.30) 


et la condition d’entrée 


u‘(.,)=ui(,) sur {(f,x)E[-1,oo[x0Q]|xeT*(r)}, (1:31) 


T£ (9 ={xe 00! v(t,x)- nx <0}; 


le choix de £ €]0,min(1, D) , de v. et de uf (voir (1.23), (1.24), (1.28)) nous garantit que 
F£(9) =T-(#) pour t z 0 et que, si u{(t, x) # 0, alors xe T£ (#). 

Nous allons considérer le problème (1.29)-(1.31) comme problème approché de 
(1.6)-(1.8). 

Comme v, est une fonction suffisamment régulière, pour tout (f,x) €] — 1,00[ xQ 
il existe l’unique caractéristique y = y‘ telle que (f, x) € y. Cette caractéristique y = 
(Yo, Ÿ), Ÿ = (Y1:::",Yn), est définie par les équations différentielles 


ee. = LT = v(#,70)) = ve(y() (1.32) 
at! 0 , at! E , € . 


avec les conditions initiales 


Yon =t, YD=x 


Si on pose 


s=inf{f <t|y(f)€e]-1,00[xQ}, (1.33) 


alors, en tenant compte que {xe 0Q|n,:v4(5,x) <0}=T£(s),onaou“s=-1ety(s)€ 


Q” ou “s>-1ety(s) e € (s)” (dans ce dernier cas on a nécessairement s > —€); le cas 
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“s = -1 et ÿ(s) e T£(—-1)” est exclu par le choix € < 1. Dans tous les cas la solution u* 
peut être exprimée par la solution de l’équation différentielle sur la caractéritique. Plus 
précisément, u° = u*(y(:)) est la solution de l’équation 

d € € € 

—U" + GU° = fe + Uj0e (D) (1.34) 


dt 


sur y = y‘ avec la condition initiale 
u° ( S) = . (1 .35) 


On remarque que, si f z 0, wo = 0, wi > 0, alors f. z20, uf z20, uŸ 20 et donc la solution 


u* de (1.34) sera elle aussi non-négative. 


Nous donnons un lemme technique. 


Lemme 1.4.1. Soient y = y‘ une caractéristique définie par (1.32) ets l'instant défini 


par (1.33). Si on pose 


À. (5) = f'euar, 
—] 


on a, pour toutt>s, 


t 
f'etaermiar <a -a élue (1.36) 
S 
QuË(s,x)1< A (s)(hs, HS Io(ts-e*,15+e1+1< 00) + (1.37) 


+(1— A(s,x) [lt ILzooqts-el+,1s+e1*1x 00)» 


où [T]* = max(r,0) (ici u1 se considère prolongée par 0 pour (t, x) € [-1,oo[ x0Q tels 


quexgT-(#)). 
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Démonstration. On a 


t t 
Jeter [octariuslie, 
S S 


ra 
d’où, compte tenu de la relation f e.(t')}dt'<1-—21;(s), on obtient (1.36). 
S 


D'autre part, on a 


max(0,5+E€) 


(CPETONEI lat, Dlee(s- dr < 


max(0,s—-€ 
< || zo(ts-e]*,(s+e1*1x 00) Àe(S); 


analoguement on a 


[GC x) + Le CONS NA Ilzo(is-e1t (sr 1x 00) Àe(S). 


De ces inégalités et de (1.28) on déduit (1.37), ce qui achève la démonstration. 


1.5 Estimation des solutions approchées sur les caracté- 
ristiques 


L'estimation de la solution approchée u‘ et de son gradient sur les caractéristiques 
y = 7°, qui joue le rôle central pour la démonstration du théorème 1.3.1, est donnée 


par le lemme suivant. 


Lemme 1.5.1. Sur chaque caractéristique y = y, la solution u* du problème approché 


(1.29)-(1.31) vérifie les inégalités 


VAOIE (acts [RO ITS Irooqrs-ert,15+e11x 00) + (1.38) 


ïf 
- J  I&G&)ldt 
+(1— h(S) Ille qs-a ts+etxam] + (À — 2e (9))| US lLeço) ent ne 
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t 
. fict)dr" 
: Î PO dt, 


max(—E,s) 


[Vu*(#)| < (1.39) 


= ae() [LACS I VHS Izoqts-et,15+e1+1< 00) + À — R(S) IV lle ts-e1+ 15+e1*1x 00) + 


fs 


. - J  o.(t)dt' 
+IVA(s)I ui — Ui ILoo(ts-e*,1s+e11x00)] 0/75) MATE ++ 


fé t 
we (t')at" 


J Î (IV) + lu (IVe (t1)e" dt’. 


max(—E,s) 


où W£(t) =. max  Oiice + ôxi Ve, jl et h(s) = h(y(s). 


1,]=1;::, 


On précise que, comme dans l'énoncé du théorème 1.3.1, dans (1.39) la compo- 
sante normale de Vu est donnée dans (1.12) (si v-n, < 0) et celle de Vu! est définie de 
manière analogue (voir (1.43) en bas pour son expression explicite). On précise aussi 
que dans les seconds membres de ces inégalités on considère u prolongée par 0 sur 


[0,c01*8Q\ U  [{r}xT-(r)] comme nous l'avons posé dans la section 1.4. 
0<1<oo 


Démonstration. Rappelons que la solution u° du problème (1.29)-(1.31) peut être construite 
par la solution du problème de Cauchy (1.34)-(1.35) sur chaque caractéristique y = y“. 


Comme l'équation (1.34) est linéaire, on obtient immédiatement, si y(s) € 0Q (et donc 


S > —E€), 
” - fictar : : : ; | fleur" ; 
lu" (n1<lu;(y(s))les + [ira +uorémeetre" dt, 
S 
et, si y(s) = y(-1) e Q (et donc s = -1), 
À , | , Jiceu")at" ; 
PO < [ira + Eye le" dt. 
= 


D’après le lemme 1.4.1, on en déduit l'inégalité (1.38). 
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Concernant la deuxième inégalité, en dérivant par rapport à x; les deux membres 
de l'équation (1.29) et en posant w = 0xU°, Le,i = Ox; fe — U*Ox; Ce et W5 ; = Ô;;,ug pour 


i=1,-.-.,n,onobtient 
n _ 
Ou + ve: VU + D (O;jCe +Ôx, Ve,j)W} = Leit+Ww6;0e(@) dans[-1,f1xQ. (1.40) 
i=1 
On considère cette équation avec la condition initiale 
w?(-1,)=0 sur Q (1.41) 
et la condition d'entrée 


w?(,)=0,u) sur [-1,f]x00. (1.42) 


On précise que dans (1.42) w° désigne la ième composante du gradient Vu composé 
par sa composante tangentielle qui est le gradient (naturel) Vr_ (y uf sur la variété T- (5) 


et la composante normale définie par 
€ 1 € € € € € 
x: Vui)|r_(p = as Ve Vr_(oUi — Oui — UV - Ve + geUï + fe + UG0e), (1.43) 
Ee°rx 


OÙ Ur,e = Ve — (Ve Nx)Nx. 
Sur chaque caractéristique y = y° l'équation (1.40) se réduit à l’équation différen- 


tielle ordinaire 


d n 
—UWi + D (Oijce + Oxi Ve, W$ = Lei+ Wet, i=1,,n. (1.44) 


dt 
W£(0 = 4) 5 we. 
i=1 


En multipliant les deux membres de (1.44) par w° et en faisant la somme de ces équa- 


ji 


On pose 


Imane BAZINE 22 Université 8 Mai 1945-Guelma 


Chapitre 1. Estimation dans W} de la solution de l’équation de transport … 


tions, on obtient 


1 d n n n 
SW) = Ÿ) (ôijce+0x, Ve,j) WjWÿ + D Leiwi+0e(D ) uw; (1.45) 
i,j=1 i=1 i=] 


cette équation différentielle ordinaire sur une caractéristique y doit être envisagée avec 


la condition initiale 
WE(9=IVu(y(DI=W,, pour se[-1,f}, 


où s est défini par (1.33) et Vu’ (y(s)) doit être entendu dans le sens précisé après (1.42). 
Si on pose 


We(f)= max loijCe+Ox, Ve, il, 
i,j=1l,.,n 
Le = (Le,1° °° »Le,n), 


de l'égalité (1.45) on déduit 
Idees Er y2 € see 2 
DIN ())° < we (WE) +1L.(DIW° (D + 0 W°(E) 2e OU . 
Donc, pour tout Ô >0,ona 
(WE(H)2<Y5(#), pourt>s, 
où Y5() est la solution du problème de Cauchy 


1 d L 
> 25 60) = (0) F5) + LD Y5(E) + QE) V Y5(t) DAUAUUDE (1.46) 
i=l 


Ys(s) = (W° +6). (1.47) 


Or, comme Y5(f) est strictement positive, pour tout f > s, on peut diviser les deux 
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membres de (1.46) par 4/ Y5(f), de sorte que l’on a 
W*(D <ys(f) pourt>s, 
où y5(f) est la solution du problème de Cauchy 


d n 
ai7o0 = We()y5(t) + ILe(D1+ 0e(0 1] > Gi ;(7(0))?, 
i=l 


ys(9 = W° +6. 


Comme y5(f) tend vers 


t t t 
__- Jœo.(r)dt! n Joet(t)dr" 
y() = We + fazer OP TTO dt’, 
rl 
S 
quand ô — 0,ona 
L A t 
= foe(t')dt! n fo(t")at" 
W(D < Wes + [azur a et / Er rene dt. 
= 
S 


D’après le lemme 1.41ona 


LÉ 
Î pur] D GS ,D}2dt < (1 Ae(5) IVuSIIze (a), 
i=1 


2. . : : 
W,< A (S(R(S IV lose ts+e]t1x00) + (A — h(S)) Vi ll Le (set 15+e+1x00) + 


+IVR(OIRÉ — lle qis-el* (s+elt1x 00) - 


t t __ 
Comme l'intégrale f : dt’ peut être remplacée par  f -dt'(voirla définition de W, 
$ max(-—E€,s) 


L:(®), 0:(P)), on en déduit l'inégalité (1.39). 
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1.6 Démonstration du théorème 1.3.1 
Démonstration. On pose 


L 
JS Gt@")+IDre (ar 
(emax(-e,s) IVc()1+1c()1+1Dve(#))dt 
\e ( t) — emax(e,s) x 


. _ 
X [Ae(S) (RS IU Il roqts-et (rer Go + À — A(SD I Îlzooqts-e1+ tre +1, 60) + 


+IVRGOIIUT — wi llre sert (s+er1xa0)) + À — 2e (DIU lle con + 


t 


L Î (fe GI+IV HN) dr]. 


max(—E€,s) 
En faisant la somme de (1.38) et (1.39) sur chaque caractéristique y = y‘ et en appli- 


quant le lemme de Gronwall, on trouve 
[uf(H)1+IVuf(n|< AP. (1.48) 


D'autre part, si on pose 5, = [s]* +e et si on considère r et r2 tels que 5: <r1 <r2<t, 


alors, en raisonnant comme dans la démonstration du lemme 5.1,ona 


r2 r2 
J'ic(ldt ’ ; 
Lu (r2)| < e Quttrte firuiar, 
11 
r2 
J'oe(')dt' Ê | ; 
IVu*(r2)l < e (vucroi+ avai etcive.(bdr) 
1 
et par conséquent 
r2 
rm JUce(@)}+IDve(r))dt! 
J(e'1 IVc()+Ice()1+1Dve()Ddt' 
[u*(r2)1+IVuf(r2)| < e1 x 


(lue + 1attrot+ farce v rca). 
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Comme cette inégalité est valable pour tout (r1, r2) € [se, £] x [Se, fé] avec r1 < r2, d’après 


le lemme 4.3 de [ASC 14],ona 


t 
JAVæeI+Ic()1+1Dve. (dt 
[uf(H)1+]Vuf (fl < ef x (1.49) 


Ë 
(lu Iat(ot+ [RIVE 


En tenant compte que l'inégalité (1.48) est vérifiée aussi pour f = 5, de (1.49) on peut 
déduire 


t 
: : SV UNI EI+IDve(#)D dt" 
[u (I +IVu (D < ee x (1.50) 


t 
x(Aë(Se) + L (AG IV GDdt"). 


Si on examine l'expression de A;(s.), on peut déduire de (1.50) qu'il existe deux 
fonctions continues v,(-) et v2(-) définies sur [0,min(1, D) telles que v1(0) = v2(0) = 0 
et que 


t 
: : Java )I+Ice@)1+IDv.(t)Ddt' 
lu (HI +IVu(DI< ef x (1.51) 


V1(E) T —E€ T 
Xe PTS RIT rose us on + À = ADI WI rose pere, Go + 


+IVR(SIITS — 2 llzeqs-et 1s+et1x00)) + A — AS) IA lys co] + 


L 
+ far +IVf- Dar) + ve). 


Compte tenu de la relation 1,(s)a + (1—1;(s))6 < max(a, 5) pour a > 0, B z 0 et des 
définitions de w}, uf, fe (INWi — willre — 0 pour € — 0 grâce à (1.18)), on déduit que le 


second membre de (1.51) est majoré par 


e"1)max(A(#), sup essB(f: 5)) + v2(e), (1.52) 


O<s<t 


LA 
: Let nb oo + Do (Isa) dt" 


t 
re (lolo + LP odt) 
0 
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A 
PUeCEg çay + Do sega) dt 


t 
B(;5)=e (aa ro + ] MPa 
S 


Cela étant, on remarque que {u*}...0 est uniformément bornée dans L®(0, 5; W1 (Q)), 
ce qui, en vertu de (1.29), implique que {ô;u*}:>0 est uniformément bornée dans 
L4(0,1,L®(Q)). Par conséquent, il existe une suite {€; 1 et une fonction 
u € L°(0,5WX(Q) n W,(0,5L%(Q)) telles que €; — 0 pour j — oo, u‘i et Vufi 
convergent faiblement-+ vers w et Vu dans L®([0, £] x Q) et ô,u‘i converge faible- 
ment vers 0,;u dans L4(0,t; L°(Q)) (dans le cas q = 1, en établissant d’abord (1.54), 


on peut procéder de manière analogue à la démonstration du point (b) du lemme 4.4 


de [ASC 14]). Donc, en faisant tendre e ; vers 0 dans (1.52), on obtient 


OIPAS < max(À(t), sup essB(f:5)). (1.53) 


O<s<t 


En appliquant l'inégalité de Hôlder aux intégrales dans l’expression de À et de Beten 
rappelant la définition de c(f), on parvient à (1.19). 


De plus, ufi vérifie, pour toute @ € €! ([-1, f] x Q) telle que 


pU,2=0 sur ([-L0]xT,(0)U( Ü UHxT,()}u (A »xO), 


O<t<t 


l'égalité 


T 
. Ô | 
Jar fu + Ve, VO +(V-ve)P— Ce;P)+ fe;P + Uÿ/ be; (DPIdx = 
—1 Q 


T 
= [ four. waras. 
-10Q 


En passant à la limite pour e; —0,ona 


t 


n 
far [rue +v-vo+gp+roiax- | Î W-mupdsat- | up, wax. (1.54) 
0  Q 0 T_(n Q 
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Comme dans (1.54) @ est arbitraire dans la classe indiquée, on peut affirmer que u 
vérifie l'équation (1.6) et les conditions (1.7)-(1.8) presque partout dans les domaines 
relatifs. 

Si ulll et ul2l sont deux solutions du problème (1.6)-(1.8), alors U = ul1 - ul] sa- 
tisfait à l'équation (1.6) avec f = 0 et avec les conditions initiale et d'entrée nulles, de 
sorte que, compte tenu que v-n >0 sur l',;(f), par le produit scalaire de l'équation (1.6) 


avec U même, ona 


1 d 


1 1 
I? =; | Ulv: as+ | = ={(V:»)) UF dx 
> qi) l 72 2 7 | CL E > v))| X 


1 
AR IOTIAET 


= 0, on obtient |U||*,.. = 0 pour tout 


» x. ., 2. 
d'où, compte tenu de la condition [U|7 T2(0) 


2(Q) | t=0 


t>0, ce qui démontre l’unicité de la solution du problème (1.6)-(1.8). 


En outre, si f > 0, uo > 0, u > 0, alors, u*i étant elles aussi non-négatives, en pas- 


sant à la limite,onau>x0. 


Remarque 1.6.1. De manière analogue à la démonstration du théorème 1.3.1, on peut 
donner une estimation de u et une de son gradient Vu dans la norme L® par les inéga- 


lités ci-dessous 


Iu(t)|roça) < max(A1(r), sup essBi(t;s)) Vtel[0,r], (1.55) 
O<s<t 

IVu(t)|ro(a) < max(A(f), sup essB(t;5)) Vte [O, £], (1.56) 
O<s<t 


t 
JA) Lo (o+I(V-v)(#)ILoo(a))dt' 
0 


LA 
Pr [lol zeeco + L LC Ireco dt", 
0 


t t 
JAI zoo (a +I(V-v) (zoo) dt" 
Bi(t;s)=es [ze CS) zocr_ (5) +f If Leo dr], 
S 


L 
J'AgE) zoo + (V0) (Leo + Dv(#) Ie ça) dt 
Ab(f) = eù (Vol 10 + 
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t 
+ fav retieo + [ue IV -V:v) (lle co)dt), 
0 


t 
J'AgE) zoo ça +l(V-v)(#)1100ço + Du (Ir (o)dt 
B2(t;s)= es UMAOI TAN: 


t 
+ [avec + [ul V(g- V:v)(1) Leo) dt). 


Corollaire 1.6.1. Soient u une fonction qui vérifie l'égalité (1.54) et u' une autre fonc- 
tion qui vérifie la même égalité où v, g, f, uo et u\ sont remplacées par v', g', f', us 
et u' respectivement avec la condition v-n, = v'-n, sur [0,t] x 0Q. Alors pour tout 
te[0,flona 


Iu(0) — u' (from) < max(A*(#), sup essB*(f;s)) (1.57) 


O<s<t 


et, dans le cas où u1 = il, on a 
Iu(e) — u’(f)Iroço) < A*(D, (1.58) 


où 


t 
JUS Lo ça +1(V-v) (Lo (a) dt’ 


A* () = [lo — upllzeçoy+ 


t 
+ [OU Piero + Ave) — v'(#)Izoço Vu (f)1e(0)+ 
0 
+(V:v-g)()-(V-v'- 8) (lie (EP) IILe(çm)dt'], 


t 
JAI Lzoo co +1(V-v) (zoo) dt’ 


BE S=e [llze1 CS) — 4 (11 cr ç9) + 


ÏÉ 
+ fauru — F6) 1o0ço) + 110) — v'(#) zoo Vu) 1(o)+ 


+08) -(V:v'- 8) (Elie col (#)ILe(çm)dt']. 
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Démonstration. Comme u — u’ satisfait à l'équation 


Ô(u—u!) 


> +V:((u-u)v)=g(u-u)+(f-f)-{(v-v)-Vu+ 


[(V:v-g)-(V:v'-g'lu, 


en appliquant (1.55) à u — u', on obtient (1.57). 


Dans le cas où 1 = u’, on voit aisément que B*(s) < A*, d’où (1.58). 


Montrons maintenant une variante du théorème 1.3.1. En utilisant les notations 
Q* = Qx]0,+oo!, V = (V, 0m) = (Ox,°*" 0x Om) et D = (V, Um) ER! xR, on considère le 


problème 


Ou+V-(uD)=gu+f dans [0,xQ*, (1.59) 
u(0, x, m) =uo(x, m) dans Q*, (1.60) 
u(£, x, m) =ui(t, x, m) sur  (J {}xT_(, (1.61) 

O<1<T 


où L_(#) = {(x, m) € 0Qx]0,00[|n,: v(t,x, m) < 0}. On suppose que les conditions A, 
B, (1.16), (1.17) et (1.18) sont vérifiées avec les modifications formelles consistant au 


remplacement de Z:, par 


Ée, = {(#, x, m) € [0,0o[ xQ x [0,001] dist((#, x, m), Ë) < €} 


avec 


= {(1,x, m) € [0,co[x0Qx ]0,00[| —€1 < v(f,X, M): nx <E1}, 


del _(#)uT (fo) par T_(t)UT, (to) (avecT , (to) défini de la manière naturelle), de la ca- 
ractéristique y par la caractéristique Ÿ définie de la manière naturelle dans [0,0o[ xQ* 
et de L°( U [{t}xT-(#)]) par L°( U [{4 x T_(9]). Pour ce problème on a l’affir- 


O<t<t O<t<t 
mation suivante. 
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Corollaire 1.6.2. Soit t > 0 (arbitraire). On suppose que 


DeL1(0,5WS(Q*;R"*D), V:5eL4(0,5:WL(Q*)), 


g, fEL TO, EWE(Q*t), ue WL(Q*), 
Umtt,x,m)=0 si mel0,molU]Mo,ool, 
uo(x,m)=0 si mel0,molU]Mo,æl 
utt,x,m)=0 si(x,mel_(t) et mel0, mIU]Mo,œl, 


HU 11) 0 si m €]0, molU]Mo,oo, 


(1.62) 


(1.63) 


(1.64) 


(1.65) 


(1.66) 


(1.67) 


où 1 < q <o et Mo > mo > 0. Alors le problème (1.59)-(1.61) admet une solution u et 


une seule dans L°(0, 1; W,(Q*)) n W,(0, 5 L°(Q*)); l'estimation (1.19) reste valable 


en remplaç ant Q, V, v par Q*, VAT respectivement et on a SUpp u(f, x,-) € [mo, Mo] 


pour tout (t, x) € [0, £] x Q. De plus, si f, uo et u1 sont non-négatives alors u l’est aussi. 


Démonstration. La condition (1.64) implique que les trajectoires qui passent par une 


partie de Q x [mo, Mol] restent toujours dans Q x [mo, Mo], ce qui nous permet de procé- 


der, dans le domaine Qx]mo—Ee, Mo +Eel (avec un € > 0 suffisamment petit), de manière 


toute analogue au raisonnement de la démonstration du théorème 1.3.1. 


D'autre part pour mé [mo, Mol], le problème (1.59)-(1.61) se réduit à 


du+V:(uv) = gu dans [0,f] x Q, 
u(0, x) = 0 dans ©, 
u(f,4)=0 sur [tx T_(@], 
O<t<t 


ce qui implique que u = 0 dans Q x (]0, molU]Mo,oo). 


On pourrait affaiblir la condition sur le support des fonctions sur [mo, Mol, mais 
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dans ce travail nous nous limitons à cette version de variante, que nous allons utiliser 
dans la section suivante. 

De manière analogue à la démonstration du corollaire 1.6.1, on peut avoir une in- 
égalité similaire à (1.58) pour deux solutions u et u”’ du problème (1.59)-(1.61) avec les 


conventions de notation analogues au cas (1.58) 


Iu(9 — u'(DIlzoças) < À*(, (1.68) 


t 
- JL Izoo (a+) +1V-D KI Loça+)dt' 
A*(t) = e 


! 
[lu — uollzo(a+)+ 


t 
+ fauru — Pie ot + 104) - (4) Ie ça Vu (Lea) + 
0 


+10:2-9(4)-(0.0 - 8) (lieu) zecs)dt|. 


1.7 Conclusion 


Il nous semble que les hypothèses du présent travail sont suffisamment faibles (au 
moins par rapport à la méthode utilisée), même si nous ne pouvons pas exclure l’éven- 
tualité d’affaiblissement des hypothèses; il nous semble que l’éventuel affaiblissement 
d'hypothèses exigera une nouvelle argumentation assez consistante, exigeant un nou- 


veau travail. 
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Application à un modèle pour 
l'atmosphère avec la transition de phase 


de l’eau 


à 7 . Es 2 . . . a . . Es 
OUS allons étudier un système d'équations quasi-linéaires, qui correspond à une 


partie du système d'équations du modèle du mouvement de l’atmosphère avec 


la transition de phase de l’eau proposé dans [SEL 11] (voir aussi [BEN 18]). 
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Chapitre 2. Application à un modèle pour l’atmosphère 


2.1 Position du problème et résultat principal 


Nous considérons notre système d'équations pour les fonctions inconnues o(f, x), 
T(t, x), o(t,x,m) dans Q* = Qx]0,0o[ où Q est un domaine borné de R° muni de la 
frontière 0Q de classe €! (1e [0, #], xe Q, m > 0); les fonctions o(#, x), x(t, x), o(t,x, m) 
représenteraient, dans le modèle de l’atmosphère, la densité de l'air sec, la densité de 
la vapeur d’eau et la densité de l’eau liquide contenue dans les gouttelettes de masse 
m respectivement. Plus précisément, en désignant par v(f, x) la vitesse de l'air (c’est-à- 
dire, vitesse pour 0 et x) et par w(f, x, m) la vitesse des gouttelettes d’eau (c’est-à-dire, 


vitesse pour ©), nous considérons le système d'équations 


d0+V:(pv)=0, (2.1) 
On +V-(rv) = -Hy](r,0), (2.2) 
00 +V-(OWw) +0m(Mmhg10) + Soft, 0)0 = S1(x,0), (2.3) 


Hgl(, 0) = =) | simo(m)dm, 
0 
het; M) = S(M)(x —- Ts), 


So(r, 0) = —hgr+ 81 (MIT Toys] + m [ Bt m'ho(m')dm', 
0 


S1(7,0) = go PCNR Ft + | Büm.m-mhotmotm- md 
0 


Tys € L4(0, f; WL(Q)), S, 81; So € WL(R), ro € WL(R: X R+), SZ 0, 80 2 0, B 2 0, tandis 
que P(:) est un opérateur lipschitzien défini sur W1 (Q*) à valeurs dans W1 (Q) tel que 
P(o) > 0 pour tout o > 0; nous supposons que 1 < q < o et qu'il y a deux constantes 


Ma, M, telles que 0 < my < M, < oo et que 


Btm',m")=0 sim +m"z>M, supp go € [Mas Mal. 
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Du point de vue du modèle du phénomène physique, 7,, représenterait la den- 
sité de la vapeur saturée (qui dépend sensiblement de la température), s(m)(7 - fus) 
la quantité de la condensation ou évaporation de la vapeur d’eau sur les gouttelettes 
d'eau de masse m (exprimée par rapport à l’unité de masse) et donc H,7(7,0) serait la 
quantité totale de la condensation ou évaporation; g1(m)[7 - 7,4] devrait représen- 
ter le taux de disparition de gouttelettes de masse m (avec m proche de m,), tandis 
que go(m)P(o) [Tr -7,$]* serait le taux de création de gouttelettes de masse m; la fonc- 
tion B(m, m') désigne le taux de coagulation d’une gouttelette de masse m et d’une de 
masse m' (pour les détails du modèle, voir [SEL 11]). 


Nous considérons le système (2.1)-(2.3) avec les conditions initiales 


0(0,)=00eWL(Q), p0o>0, (2.4) 
x(0,)=70€WL(Q), 702>0, (2.5) 
o(0,)=00€eWi(Q*), 0020, suppool(x,) € [Ma Mal VxE Q, (2.6) 
et les conditions d’entrée 
Olrw(n = @1» @12Z0, (2.7) 
Hlrw(p = T1, 1120, (2.8) 
OlFw(p = 01 0120, (2.9) 


suppo1(f,:,:) Le 0Q X [Mas Mal VtE [O, #1], 


où F9) = {x € 0Q ns: v(E,x) < 0}, FM) (9) = {(x, m) € 8Qx]0,00[| 74: w(t, x, m) < 0}. 
On suppose que les conditions A, B, (1.16), (1.17) et (1.18) sont vérifiées pour w1 = pi et 
ui = 7, en remplaçant L_(#) par LV) (#) et pour w1 = 01 en remplaçant v par wet T_(f) 


par ['(#)(#) (et avec des modifications formelles évidentes comme dans le corollaire 
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1.6.2). On suppose aussi que 
veLT(0,WI(QR)),  V-veL1(0,: W1(Q)), (2.10) 


we Ld(0,EWL(QY;R),  V:weL1(0,5 W1(Q*). (2.11) 

Le résultat principal de cette partie est le théorème suivant. 
Théorème 2.1.1. Sous les hypothèses posées ci-dessus, il existe un t* €]0, t] tel que le 
système d'équations (2.1)—(2.3) avec les conditions initiales (2.4)-(2.6) et les conditions 


d'entrée (2.7)-(2.9) admette dans l'intervalle de temps [0, t*] une solution (e,x,0o) et 


une seule telle que 
pe, re L(0, 6"; We (Q)) n W, (0, £";L°(0)), 


oeL*(0, 1"; We. (Q*)nW;(@,r";L°(0*)); 


etonapz0,r>=0,02>0et 


suppo(f,x,)C[Ma Mal V(t,x)E [0,1] x Q. 


2.2 Étude du système des équations linéarisées 


L'équation (2.1) étant linéaire, l'existence et l’unicité de la solution résultent du 
lemme 2.2.1 montré ci-dessous. En ce qui concerne les équations non linéaires (2.2) 
et (2.3), on examine d’abord les équations linéarisées et puis on cherche un point fixe 


d’un opérateur défini par la solution de ces équations. 
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On considère deux fonctions données 7 et © et le système des équations 


d0+V:(pv)=0, (2.12) 
On+V-(rv) = -Hy](r,0), (2.13) 
00 +V-(OWw) +0m(Mmhe(T)0) + So(T,0)0 = S1(7,0), (2.14) 


qui est le système linéarisé de (2.1)-(2.3). 


Lemme 2.2.1. Six € L°(0, 1; W1 (Q)) er o € L®(0,1;, W1(Q*)), alors il existe une so- 
lution (0, 7x, o) et une seule du système (2.12)-(2.14) avec les conditions initiales (2.4)- 


(2.6) et les conditions d'entrée (2.7)-(2.9), solution dans la classe 
e,xe L°(0,5WS(Q)nW,(0,5L%(Q)), 


oeL*(0,5W.(Q*)nW,(@,5L°(Q*); 


de plus, sir eto sont non négatives, alors 0, 7 et o sont aussi non négatives. En outre, 


il existe une constante C > 0 indépendante de v, w telle que les inégalités suivantes 


soient vérifiées 
IQ ro t0#wz, con < MAX (A0, sup essB,(s)), 2/15) 
O<s<t 
(TI zoo rw, çan < MAX (A7, sup essB;(s)), (2.16) 
O<s<t 
O1 (0,7 w2, (a+) < max(4,, sup essB(s)), (2.17) 
O<s<T 
où 


qg-1 
LE 4 (DvlrawrrecoytlV-2l;a0 rw 
4, = € ( LA (0,5 L%(Q)) eo on) IlColl we, cop 


q-1 
Œ-s) T (IDvlzaçsrrooçon tV-vl;a 0 rw 
B,(s)=e (IDvlzaçrrcocy co) 11) croco 


q-1 
5 4 (IDulzsosrom+lV-vl;«0zw1 cor) 
Ar = € L4(0,FL%(Q)) 14(0,:wL(Q)) x 


= 
q 


Ct|o!| =wl _ _ 
Xe POELE SE [lToll WL (Q) +C I TvsllLa(0,7:w2, 0» 10 10 r:w2,(0+))1> 
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eo Ce = : 
Br(s) _ Un q (lDolzacarsecy+IV-vl1a rm co) +CE-9 1O1 00 (5 7:w2,(0+) " 
_ gi = 
x [||71(5) Île ra cs) HOUR S I Tvsllza(s rw, co» 10 Loos, rw, co+l> 
LI de 
= PL (Du saozze y +17 w leu cry vs 0 wa on) ; 


CT(IT = +15 = ie 
re lornnan tel araen) one qe + CET us ow et 


Flo Eme co) + Ole meet) +0le rw co») 


= 
LS a = : # = : 
Be (S) = _ S) (Dwlbss ee, IF ul Lacan cart + ClvsIlLa tv to) | 


re ecrnten leusbiat) for (sy us + 
ne 


= ga = = 
+C((Œ- Se I TvsllLa(s rw, co) HL= STI zo0(s, rw, 0) + 
+(Œ—S) IT zoo(s, rw, +) A + IC le(s Ewz, c)))- 


De plus, si on a 
suppo(f,x,:) [ma Mal pour tout (t,x) € [0,f] x Q 


alors 


suppo(t,x,:) [Ma Mal pour tout (t, x) € [0, £] x Q. 


Démonstration. On considère l'équation (2.1) avec la condition initiale (2.4) et la condi- 


tion d’entrée (2.7); afin d'appliquer le théorème 1.3.1 on prend 
U=Q, Uo=@o W=Er, f=0, g=0; 


d’après l'inégalité (1.19) on obtient l'inégalité (2.15). 


On considère maintenant l'équation (2.2) avec la condition initiale (2.5) et la condi- 


Imane BAZINE 38 Université 8 Mai 1945-Guelma 


Chapitre 2. Application à un modèle pour l'atmosphère 


tion d’entrée (2.8); pour appliquer le théorème 1.3.1 on prend 


U—T, Uo—Ho, Ui — 71, 


CO [ee] 


g= - [smaumam [= us | sm uma. 
0 0 


En vertu des conditions posées sur 7,, et s(m), on voit qu'il existe une constante posi- 


tive c, telle que 


18 lu, co £ C11lO lu, co) 


et 


FAT = CT vsllw co 10 lu, co): 


En utilisant l'inégalité de Hôlder, on a, pour 0=< 5 < f, 
= 421 _ _ 
ES) 1 1Slzas rw con < CE HIT Le(sFWL (a+) 


et 


D ai Lai L 
ES) 1 Mflrasrmgoy € CES 4 IToslra(sews oO lLe(s rw a+) 


avec une constante C. Donc, d’après l'inégalité (1.19), on obtient (2.16). 
On considère maintenant l'équation (2.3) avec la condition initiale (2.6) et la condi- 


tion d’entrée (2.9) et on prend 
U=O, U9=00 U=O0p V=W, Um=Mhgl(T), 


g=—So(r, 0), f = S1(7,0). 


En vertu des conditions posées sur 7»s, S(M), go, g1, B, P, on constate qu’il existe une 


constante positive c: telle que 


Sr ço+ < C2(1T 1 y Q + | Tosll y Q +10 ll y o+) 
(A+) eo (Q) eo (Q) (A+) 
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1 Flu co+) = C2(IT we co Le IT vs1l we, co) + 1012, ça+)) Lu 101 ça+))- 


En utilisant l'inégalité de Hôlder, on obtient 
D ai - _ _ 
Es) 7 IglrsrwLeoty < CE - UTILE we con + MTL (sm (0) + 


_ AT, -— 
+(t—s) 4 ITvsllza(s rw, 0) 


et 
= a=1 _ g21 = _ 
(84 IF zacszwz,o+) = C(Œ- s) I Tvsll Lacs rw, co) +E-$ IT Loo(s, rw, 0) + 


+ DIT rmg en) + 10e sn tot) 


avec une constante C. D’après le corollaire 1.6.2, on en déduit (2.17). 
Pour la dernière assertion du lemme, on remarque que, si suppo(f, x, +) € [Ma, Mal 


pour tout (f, x) € [0, é] x Q, alors supp f(f,x,:) € [Ma, Mal pour tout (f, x) € [0, f] x Q, 


donc le résultat découle immédiatement en appliquant le corollaire 1.6.2. 


Pour utiliser le lemme 2.2.1 pour la démonstration du théorème 2.1.1, nous intro- 


duisons 
O<s<t 
A0) = 2max(||0ol w2 (a+ Sup essllo1llyz Œurçp)) + 1: (2.19) 
O<s<t 


Alors on a le lemme suivant. 


Lemme 2.2.2. Il existe un ti €]0, t] tel que, sir et o appartiennent à L°(0, ñ; W1 (Q)) 


et L°(0, ñ; W1(Q*)) respectivement, sont non-négatives et vérifient les inégalités 


[T1 Le (0,1 :W21,(0») = Ar] (ti), 110,1: (0+) = À 0] (ti), (2.20) 


alors la solution (0,7,0) du système (2.12)-(2.14) avec les conditions initiales (2.4)- 
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(2.6) ef les conditions d'entrée (2.7)-(2.9) vérifie les inégalités 


[T1 (0, n:W2,(0)) £< At), 101110, 1:w2,0+) < Ajo](ti) (2,21) 


eton «a 


suppo(£, X,+) [Mas Mal V(E,x) € [0, ti] x Q. (2.22) 


Démonstration. On voit que, si on choisit r suffisamment petit, dans (2.16) et (2.17) le 
deuxième membre peut être inférieur à Arr) et à Aj(t) respectivement, ce qui nous 


permet de déduire (2.21) sous la condition (2.20). La relation (2.22) découle immédia- 


tement du lemme 2.2.1. 


2.3 Démonstration du théorème 2.1.1 


Soit t €]0, 1], où f est donné dans le lemme 2.2.2. Nous définissons le sous-ensemble 
F, de L®(10, t] x Q) x L®(10, #] x Q*) par la relation : (x, o) € F; si et seulement si (7,0) 


vérifie les conditions 


ne L*(0,t: WL(Q)), oe L°(0,1, WL(Q*)), 
T,0 >0, 

suppo(#’,x,-) € [ma, M] V(’, x) € [0, r] x Q, 
TI (0, rw, (07 < Ati (6), 


lo ze(0, sw ay < Ato1() 


On considère, pour tout f €]0, ñ], l'opérateur G; : F; — L®°(10, #] x Q) x L®(10, x Q*) tel 
que G;(7,0) = (7,0), où (7,0) est la solution du système (2.13)-(2.14) avec les condi- 
tions initiales (2.5)-(2.6) et les conditions d’entrée (2.8)-(2.9). 


En vertu du lemme 2.22 on a 


G(FYH CF; Vtel0, fl. 
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Donc, pour démontrer le théorème 2.1.1, il suffit de démontrer qu’il existe un f*, 0 < 
t“ < t:, tel que l'opérateur G;: restreint à F,+ soit une contraction. 


Pour i € {1,2} on considère (7,0 %)e F, et (7,0) = G,(m,5°)). Nous posons 
P 


II _ 72 …: 70), > _ 52 - où), 


n=rx2-r0, »z=02-0û 


Comme =; =n% sur LJ [{tx1(5], en appliquant l'inégalité (1.58) à 


0<f<t 
u=7®, u =7n0), v= 7", Uo = Up = To, U = WU, =, 
2e OO 
… —(2) pa . 
g=-|stm)0" (mdm,  g'=-| s(mo (m)dm, 
L, 0 
9 OO 
Î= | s(mo(mdm,  f= Rs | s(m)o dm, 
8 0 


on obtient, en vertu du lemme 2.2.2 et de la définition de F,;, 


IT zæGo,r1xo) < An, (2.23) 


q-1 


An — Ce! q IV-v1l La (0,1:100ça7 + CEAol (D Le (2.24) 


= g=1 _ 
X 12 ze go,#1x0+) Ar) + € 4 TTusllzaco rw, cpl 


D'autre part, pour estimer ?, utilisons les notations 


Cs= sup (m+m')B(m,m'),  C;= Isle qu» 
m,m'>0 


Cu =supg(m), Cy = SUP 8o(mM). 


m>0 m>0 
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Alors, comme 00 =01=0% sur |} [{4 x T((5], en appliquant l'inégalité (1.68) à 
O<t<t 


U=0, U=0 y, W=d=0p M=W=0y 


8=S@°,07), g'=9@°,00), f=8@%°,07), f'=S(@,00), 
D=(wmham?), D'=(w,mhgm 0), 


on obtient 


1ZlIzogo,r1xat) < A5; (2.25) 


Éms = — gl = 
e! 7 Ilan) (@+MaCs+Cg Am O+E 7 Muslla( + CptAlo( 


Ax = ,EWeo() (2.26) 


— gi 
x {CglCP1ZILeGo,11x0+) (At (©) +194 UTvusllza(o, rw, op) + 
+tCp(Ao]() + DIT L+00,1x0)1+ 
+t[20 + MC; + Ca) IL g0,11<0) +2MaCgllËllego,1x0+)Ato(Ô}, 


CP étant une constante due à l'opérateur lipschitzien P(.). 

L'estimation de II et celle de Z étant établies, on va raisonner séparément pour le 
cas g > letle cas g=1. 

Dans le cas q > 1, des relations (2.23)-(2.26) on déduit qu'il existe une constante Ci 
telle que 


TT ze go,#1x0) + I Zlzego,r1xa+) < 
cu q=1 _ ” 
< Ce +0 a + (IIIe go,rtxe + IEllLegorxot)  VtE 0, ñ]. 


Donc, si on choisit t* e]0, f.] de telle sorte que 


RER a=1 
GG TRE ed 
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alors G;: sera une contraction. 
D'autre part, si g = 1, on remarque que, quel que soit € > 0, il existe un ff > Otel 


que 


IV: VÎlL1(0,1;1(0») < E; UTvsll z1(0,1::w2, co» < €, IV: W1 L1(0,1*;L°(Q+)) < E€. 


Cela étant, si on examine les expressions (2.24) et (2.26) de An et de A;, on constate 


qu'on peut choisir un f* > 0 suffisamment petit de telle sorte que l'inégalité 


(TTL 1° go,r1x0 + IZlzogo,rtxe+) < K(IHTILzoego,#1x0 + Ze Go,r1x0+)) 


soit vérifiée avec un k < 1. Donc G,: sera une contraction. 

Ainsi on a montré que dans tous les deux cas g > let q = 1; il existe un t* >0Otel 
que G,- soit une contraction. Donc il existe le point fixe (7,0) € F;+ (unique dans F;+) 
pour l'opérateur G;:. En rappelant que la solution 0 de l'équation (2.12) est donnée 
dans le lemme 2.2.1, on constate que (0,7, o) ainsi obtenu est une solution du système 
(2.1)-(2.3) avec les conditions mentionnées en haut. 

La non-négativité de p, de x et de o résulte des conditions s > 0, go > 0, P(o) >0 


pour tout o > 0, comme dans le théorème 1.3.1 et le corollaire 1.6.2. O 
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Affaiblissement de la condition pour 
l'estimation dans W} de la solution de 


l'équation de transport avec l'entrée 


3.1 Introduction 


4] DANS [ASC 14], les auteurs ont obtenu entre autres l'estimation dans wi de 


C2) la solution de l’équation linéaire, en supposant que sur la frontière du do- 
maine la composante normale de la vitesse de transport s’annule. En suivant l’idée de 


cet article, mais en admettant l’entrée et la sortie à travers la frontière, dans [BAZ 18] 


45 
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(voir chapitre 1), on a établi une estimation dans W1 de la solution de l'équation de 
transport linéaire. Or, la condition sous laquelle on a démontré cette estimation était 
quelque peu de nature technique. En effet, en considérant l'équation 


ou 
— +V: = + 
a (uv) = gu+f 


avec un champ de vitesse v (get f sont des fonctions données) dans un domaine borné 
QCR’ et la condition d’entrée 


Ut, x) = u1(E, X) 


pour (f,x) ER; x 0Q tels que n,:v(f, x) < 0 (n, désigne la normale extérieure à 0Q), on 
a supposé l'existence d’un £2-voisinage de la zone {—€1 < n(x):v(f,x) < €1} dans lequel 
les fonctions données sont suffisamment régulières pour avoir la solution classique. 

Dans ce chapitre, nous proposons d’affaiblir cette condition posée dans [BAZ 18] 
pour obtenir un résultat analogue sur l'estimation dans W1 de la solution de l’équa- 
tion de transport linéaire avec l'entrée. Pour ce faire, nous utilisons une régularisation 
particulière des fonctions données et estimons la mesure de la zone où la solution du 
problème approché peut avoir de mauvais comportements. 

Ilest bon de rappeler que durant ces dernières années, l'équation de transport a été 
beaucoup étudiée, en particulier sur l’affaiblissement des conditions sous lesquelles 
on démontre l'existence d’une solution (voir par exemple [AMB 04], nous citons aussi 
[DrP 89] pour son importance). Même si les travaux traintant la condition d’entrée ne 
sont pas nombreux, le problème a été étudié dans [BAR 70], [CRI 14], [BOY 05]. L'intérêt 
pour l'estimation de la solution dans WL est dû en particulier aux applications aux 


équations qui ne sont pas linéaires, comme il a été traité dans [ASC 14] et [BAZ 18]. 
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3.2 Position du problème et résultat principal 


Soit Q un ouvert borné de R’. On suppose que sa frontière 0Q est de classe €}, Q 


étant d’un seul côté de 00, et qu’il existe un €0 > 0 tel que, pour tout x € 00, on ait 
x-en;seQ VeEel]0,€ol, (3.1) 


où ñ, désigne la normale extérieure unitaire à la frontière 0Q au point x € 00. 
Sur Q on considère le champ de vitesse v = v(f, x) définissant le flux de transport. 


On suppose que 


ve L}, (0,00; W(G;R°")), (3.2) 
V-veL}, (0,00; We (Q)). (3.3) 
On définit 

T_(D ={xe0Q]nxv(t,x) <0}, (3.4) 
lo(f ={xE0Q ns. v(t,x) =0}, (3.5) 
T4(D ={xE0Q]nx-v(t,x) > 0}. (3.6) 

On pose aussi 
Ao = Ü (8 x To(n). (3.7) 

tz0 


Avec v et Q introduits ci-dessus, on considère l'équation 


Ô 
_ +V:(uv)=gu+f dans Q, (3.8) 


où get f sont deux fonctions données à valeurs réelles. L'équation (3.8) est envisagée 
avec la condition initiale 


u(0, x) = up(x) pour xeQ (3.9) 
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et la condition d’entrée 


u(t,x)=u(tx pour(xe UJ [{Hxr_(|], (3.10) 


0<t<oo 


uo(x) et u1(f, x) étant des fonctions données. On suppose la condition de compatibilité 


Uo(xX) = u1(0,x) pour xel-(0). (3.11) 


Pour préciser les conditions sur les données, nous introduisons les notations Vr_(nu1 
et n;Vu:le symbole Vr_(#u1 désigne le gradient de u, sur la variété l_(f) dans le sens 


naturel, tandis que n,:Vu, est défini par 


nx° VU] = (-v5-Vroti-0u-WV-v+gu+f) (3.12) 


Vehx 


avec v; =v-(v-n;)nzx Sur L_(f). 
Désignons aussi par y(1,2 = {(Yo(t),Y1(6), +: ,Yn(0))}re1,,x la Caractéristique passant 


par le point (f, x) € [0,co![ xQ se définit par le système d'équations différentielles 
2 COHEN d 1) = v;(y(E = 

arroU)=L SU), ÉCRIRE 

avec la condition initiale 


Yo =t, YjÜ=x;, j=1l,..,n; 


OÙ lmax est l'intervalle maximal d'existence de la solution. 
Pour énoncer le résultat principal, nous formulons les conditions suivantes : 


CONDITION A. On pose 


V5(Ao) = {(f,x) € [—1,00[ xQ[dist((r, x), Ao) < Ô}, (3.13) 
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E(V5 (Ao)) = {(£, x) € [0,00[*Q 1 7,2 N V5 (Ao) À GI. 


Alors pour tout t>0ona 


mes(E(V5s (Ao)) — 0 pour ô — 0. 


CONDITION B. Si (fo, Xo) € Ao, alors 
(Lo, Xo) (a Ÿ(,x) VE X) € [0,co! x 0. 


On a le théorème suivant 


Théorème 3.2.1. On suppose que les conditions À et B sont remplies et que 


geLp,. Oo We, EL, Oo We(D), ue e Wé(O), 


loc 


v, Vr_cUr fu, g—V-veL®( U[J [{xT_(9]) vVr>0, 


O<t<t 


HV LS NBI nl) 0, 


O<t<t 


uiest uniformément continue sur [tt x T_(f] pour toutt >0et 
O<t<t 


y admet presque partout la dérivée par rapport à t. 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


Alors il existe une fonction u € L® (0,00; W1 (Q)) et une seule qui satisfait à l'équation 


loc 


(3.8) et aux conditions (3.9)-(3.11) et, pour tout t > 0 on a 


Ur 7w, ça) < Max(A, sup ess B(s)), 
O<s<t 


où 
Si 
q ; : : : 
AD Uallza one cop * IV VlLa 0 we on *DVl1a0£1 ça), 


GE 
q 


x [|| Uollw1, (0) +1 MAPTATAONE 


(3.21) 
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g=1 
q 


ni UBlza sem y larme on DVI La(sE Leo») U 


B(s) = 


Gil 
q 


* [4 (Su crc + ÉD T Flzacs eme cpl: 


3.3 Formulation du problème approché 


Pour démontrer le théorème 3.2.1, nous démontrons d’abord le résultat sous une 
condition un peu plus restrictive, plus précisément, au lieu de (3.2), nous allons d’abord 


procéder sous la condition 


OO 
veL;. 


(0,00; W(Q;R°")), (3.22) 
Nous écrivons l'équation (3.8) dans la forme 


dçu+v-Vu+cu=f, c=V:v-g. (3.23) 


Nous allons considérer l'équation approchée de (3.23), en définissant l’approxima- 
tion des fonctions v, c, f, uo, 1. Pour ce faire, on commence par prolonger les fonctions 


v,C, f pour t <0 par 
v(,:) =0, c(t,)=c(-t,:), FH) =0. 


On prolonge également uw, par 0 sur [—-1,0[ x0Q et aussi sur l}(#) U lo(é) pour chaque 
t>0 (ce dernier prolongement sert seulement pour simplifier la notation). 
Pour la régularisation de fonctions, nous introduisons d’abord une famille de fonc- 


tions {h(e; t, X)}ocee1 de classe €!([-1,00o[ xQ) satisfaisant aux conditions 


h(£;t,x)=e si dist((£, x), Ao) 23e, (3.24) 
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1 1. se 
4 dist((E, x), Ao) < h(E; f, x) = > dist((#, x), Ao) si dist((f,x), Ao) < 3€. (3.25) 


Soient d](s) et é1(y) deux fonctions régularisantes canoniques ayant les supports {|s| < 


1}CRet{lyl<1} CR" respectivement. Pour chaque € > 0 on définit 


s! 


l 
el 1 = ÿ; (— ——— ; 3.26 
(xt ) h(e:t, x) 1 EL 520) 


y! 
El G)= SG) (3.27) 
î, - 
[£,x] RAGE) dx" 
Posons 
vil(£, x) =[[ DU pots, Pdyds. (3.28) 


On remarque que le choix de h (voir (3.24)-(3.25)) et la définition de GE é cn vlel 
(voir (3.26), (3.27), (3.28)) laissent invariants les ensembles T_(#f), lot) et l.(f), c’est- 
à-dire, si on pose 


TEl(n ={xedQ]n,-vEl(r,x) <0}, 


(HD ={xe0Qnx:vll(r,x) = 0}, 


r'l(n ={xe0Q ns: vll(r,x) > 0}, 


alorsona 
réO=rit IFOEmE, LIÉE: 0 Vie 
Posons aussi 
c'(£, x) = NETTS Je (x— yc(s, y)dyds, (3.29) 


fl, x) = [LL gl, (£— Ja (x— y) f(s, y)dyds. (3.30) 
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En ce qui concerne la donnée initiale wo, nous posons 
dy (2) = 00.6) [ EXP UP) dy. (3.31) 


D'autre part, pour définir la donnée d'entrée régularisée, nous définissons 


00) __ HG Eh (3.32) 
éatr,2] D : 
d Lo Eds 
où dS est l’élément d'intégration sur 0Q. Avec ceci on définit 
(2 = [ Î AE — Eau (s, NAS, ds. (3.33) 


Maintenant nous pouvons formuler le problème approché : trouver la fonction 
u(f, x) qui vérifie l'équation 
El. vu+cllu= fEl4 ue dans [-1,00[xQ, (3.34) 


ôt 


avec la condition initiale 


u(—-1,:)=0 dans Q (3.35) 


et la condition d’entrée 


u=ull  sur{(t,x)e[-1,0[x0Q|xeT*(r)}. (3.36) 


3.4 Solution du problème approché et ses estimations 


Rappelons d’abord l'existence et l’unicité de la solution du problème approché. 


Lemme 3.4.1. Le problème (3.34)-(3.36) admet une solution u et une seule. 


Démonstration. Comme la fonction vl‘l est régulière, en considérant les caractéris- 
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tiques yl1(9) = yl(to, xo; t) définies par l'équation 


t 


7 = xo + Î vll(s,yl(s))ds, (3.37) 


lo 
avec 


Goo) € (-1 x QU{ U Anxr(»))=5!, 


tz—1] 


on peut obtenir la solution u par la résolution de l’équation 
d 
ET CO) = VE TE Duty ED) (3.38) 


+gll(s,yEl(u(s, vo) + fl, El), 


avec la condition initiale 


[e] si (t7*/(H)Ee{-1xQ 
U(o,Y * (fo) = | 
ulel (to, 7 (o)) si (to, Y €} (to) eUs tx rl(6) 


Pour chaque f > 0 on définit 


Are LA X OI Ù Gtxr))n Vito) # 6}, (3.39) 


ô 
O<E<7 


où G(t,x; y!) est la restriction à [—1,f] x Q de la caractéristique y!‘ (to, xo;-) telle que 
VE (to, Xo; t) = x (avec un certain (fo, xo) € Z*). 
Rappelons une conséquence immédiate de la définition de 4,7. Pour cela, dési- 


gnons par fe , la caractéristique définie par vlEl et passant par (f, x). 


Lemme 3.4.2. Si(t, x) € ([-1, 1] x Q)\A,7 et si0 < € < è, alors on a 


h&t,1)=e Né,x)ey., (3.40) 
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où h(e; t', x’) est la fonction introduite dans (3.24)-(3.25). 


Démonstration. I résulte immédiatement de la définition de 4,7. 


On va estimer la solution u du problème (3.34)-(3.36) et ses dérivées par rapport à 
X= (x%,::,Xh) sur les caractéristiques. Pour cela, nous suivons le même raisonnement 
que le lemme 5.1 de [BAZ 18], et nous limitons toutefois aux estimations sur ([—1, f] x 
Q)\4; 7. Cette limitation et le lemme 3.4.2 nous permettent d'épargner l'estimation des 
valeurs de |u]| et de |Vul|influencées par la régularisation particulière dans le voisinage 
de A, (voir (3.25), (3.26), (3.27)). 


Comme dans [BAZ 18], définissons d’abord la fonction auxilière 


[£,x 


À (s) = Î 9! .(W)dr,  dist((f,x), A0) 23e. (3.41) 
-1 


Lemme 3.4.3. Soit y = y} une caractéristique définie par (3.37) avec le point de dé- 


part (to, Xo) € Z*. Alors on a, pour toutt > to, 


Ë 
Î YA DIdE £ (1 eo) luollzeco, (3.42) 
lo 
lu Go, 1 < Ae(to) ua ze -er ttoe" 1x 80). (3.43) 


où [T]* = max(r,0) (ici u se considère prolongée par 0 pour (t, x) € [-1,oo[ x0Q tels 


quexgT-_(#)). 


Démonstration. On le démontre de la même manière que le lemme 4.1 de [BAZ 18]. 
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Lemme 3.4.4. Sur chaque caractéristique y = y} ayant le point de départ (to, Xo) € Z* 


et telle que | 140 45,7 = ®, la solution u du problème approché (3.34)—(3.36) vérifie, 


pour t € [to, t] avec (f,y(®)) € y, les inégalités 


[u(y(O)I < OI U1|l ze (IIt-e1*,1to+e1*1x 00) + (3.44) 
t t t 
CO fil" pidr" 
+ Je" + [eue at, 
max(—E, to) 
[Vu(y(D)|=< (GT liege orerxa0+ (3.45) 
f 
J  o.{yU')dr 
+(1 — 1: (to) IVuol 12e Je" 4 
t t 
[e] ! ! [e] ! Ie ! 
: il (VF œGDI+ lu DIV E (y) e" dt. 
max(—E€, to) 
où 
we, max 60 +00 h  IVaalie=Vr-otllze + Inx-Vuilie 


(pour cette dernière notation, rappeler (3.12)). 


Démonstration. On le démontre de la même manière que le lemme 5.1 de [BAZ 18]. 
(Par rapport au lemme 5.1 de [BAZ 18], ce présent lemme est plus simple. Car ici, nous 
considérons seulement les caractéristiques contenues dans ([0,£] x Q)\4;; de sorte 


que l’opérateur régularisant reste invariant (voir (3.40)) et donc les termes dûs à la va- 


riation du noyau de l'opérateur régularisant disparaît.) 
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3.5 Estimation sous la condition de la régularité de v de 
type L® 


Dans cette section, nous allons établir l'estimation de la solution uw dans 
L®(0, f;, W1 (Q)) sous l'hypothèse de l'appartenance de v à L® (0,00; W1 (Q;R”)) (condi- 


loc 


tion (3.22)). 


Lemme 3.5.1. Désignons par ul la solution u du problème approché (3.34)-3.36) 
avec € €]0,1]. Soit t > 0. Pour chaque 6 > 0 fixé, les solutions approchées ul: 


convergent vers une fonction u dans ([0, f] x Q)\A;; pour € — 0 et u satisfait à l'in- 


égalité 
sup ess||u|| 1 (p, ) < max(A, sup essB(s)), (3.46) 

O<s<t | O<s<T 
où 

_41 

A=e ” Vlzorwgert lacune lDvlie pret) 
_4=1 
x [|| Uollw1,(o) an APT TAONE 
en = ” 
TOO 
B(s) = oi q US zacszme ot IV-VlLa (rw ont PvlLet(s, GE co), ; 
= GE 
x [|] 2 (S)Ilw cr cs) AS) RG APTE TA NE 
dans cette expression 


Ds, ={xEe Q1(,x) g A5: 
1241 CS) Il, = 
= [| u1(s)zæcr_ (5) + IVr_ (9 41(S)zæecr_ (5) + x W1(S)ze(r_ (5). 


: -1 - ie 
Dans le cas où q =, on remplace . par 1. De plus, si f, uo et ui sontnon-négatives, 


u est aussi non-négative. 


Démonstration. Pour démontrer (3.46), nous procédons par les estimations de [ul®\(#, x)| 
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et de |[Vulfl(#, x)| sur les caratéristiques y = y!‘ définies par vlfl. Or, comme toutes 
les caractéristiques y!‘ avec 0 < € < ÿ passant par Ds, ne rencontrent pas À; (voir 
(3.13)-(3.39), voir aussi le lemme 3.42), on prend en considération l'inégalité (3.44). 
Cela étant, de la même manière que dans la démontration du théorème 3.1 de [BAZ 18] 
(plus précisément, la déduction de (6.4) à partir de (5.1)-(5.2) dans [BAZ 18]) on déduit 
de (3.44) l'inégalité 


t 
[e] {el Savoir E»)+1Dv (y »Dar 
[u* (y) + [Vu (y(D)l< es : (3.47) 


X (e10) [Ae(s) I 1] go(tts-e1+,(5+e1+ 1: (600) +(1—2;(5))| Uollwz co] + 


t 
+ j fo )+IV fo )Ddt) +v2te), 


Se 
où v.(-) et v2(-) sont deux fonctions continues définies sur [0,1], telles que v1(0) = 
V2(0) = 0. 
Compte tenu de la relation 1. (s)a+(1-—1.(s)) 6 < max(a, B) pour a >= 0, B>0 et des 


définitions de vf, clél, fl£l, on déduit que le second membre de (3.47) est majoré par 


e"1)max(A(#), sup essB(f: 5) + v2(e), (3.48) 


O<s<t 


L 
ap = ot tDP sand 


t 
= (lolo + JA dt) 
0 


t t 

_ JU 2 co + D (ILæ (ay) dt’ 

B(t;s)=es 1 (ze lues cr co +f RAP RE 
$ 


Cela étant, on remarque que {ul}, 0 est uniformément bornée dans L®(0, f; wi (Q)), 
ce qui, en vertu de (3.34), implique que {6;ulfl},,9 est uniformément bornée dans 
L4(0,f;L®(Q)). Par conséquent, il existe une suite {e fn et une fonction 
u € L°(0, 1: WE (Q)) n WJ(0, 5; L°(Q)) telles que £; — 0 pour j — co, ulfil et Vulei! 


convergent faiblement-+ vers u et Vu dans L®(([0,£] x Q)\A,7) et 0sulil converge 
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faiblement vers ô;u dans L4(0,f; L°(Ds:)) (dans le cas q = 1, en établissant d’abord 
(1.54), on peut procéder de manière analogue à la démonstration du point (b) du lemme 


4.4 de [ASC 14]). Donc, en faisant tendre € ; vers 0 dans (3.48), on obtient 


ul ço < max(ÂE®, sup essB(r; 5). (3.49) 


O<s<t 


En appliquant l'inégalité de Hülder aux intégrales dans l'expression de À et de Beten 
rappelant la définition de c(f), on parvient à (3.46). 


De plus, lil vérifie, pour toute y e €!(([-1, ] x Q)\A, >) telle que 


pU,x)=0 sur ô(([-1,4 x Q)\ A3) \ Ü x T_(D)), 


O<i<t 


l'égalité 


n 
Ô 
Jar fu SE+ vlEil Vo +(V: vEil)o = clip) + fil + ul (Dpldx = 
—1 Q 


T 
= [fo nu ot mdrds. 
-10Q 
En passant à la limite pour e; — 0,ona 
T a T 
Jar fu +vvo+ go + foax= | Î W-mupdsat- | up(0, max (3.50) 
Q 0 Q 


0 [_(® 


Comme dans (1.54), @ est arbitraire dans la classe indiquée, on peut affirmer que u vé- 


rifie l'équation (3.8) et les conditions (3.9)-(3.10) presque partout dans ([0, £] x Q)\A, 7. 


Maintenant nous voulons montrer que la mesure de A; tend vers 0 quand à tend 
vers 0. Pour ce faire, nous allons d’abord estimer la différence entre yl£l(to, xo; t) et 


Yo; Xo: t). 
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Lemme 3.5.2. On a 
t 
y (to, xo; D — Yo, Xo; 01 < (2+suplvl)eL(t- +1 | (s—- te CS gs)+R(E), (3.51) 
lo 


où R(E) tend vers 0 lorsque £ — 0. 


Avant de démontrer ce lemme, en désignant par ve les caractéristiques définies 


€] 


par UE c'est-à-dire, 


t 
[e] [e] 


[e] — 
Ya (Lo, X0; D) — X0 +[ Vi} (S, Ya (s))ds, 


S+E 
[el [el - 1 ,,[€l [el UV 
er = fl V(S Ya), ts —S)ds, 


et on donne l’estimation suivante 


Lemme 3.5.3. On a 


Es 
I Go Xos D — Yo, Xo; D < ELsuplvi(£- t0+L 1 (ste Vds) FR(E) (3.52) 
lo 


Démonstration. On a 


t 
17 os Xoi D — Yo; Xo; 1 = | [. WEST) - v(s,7(D)ds| 
lo 


t t 
< | Î GG) - vs, 7 SD) ds +| Î GS, 79) — vs, 79) dsl 
to lo 


t t 
< | Î GG) - vs, 7h (S))ds|+ L L AOC 
to to 
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Maintenant, examinons la différence intégrale 


1 (ol (5, (s)) En U(s, je (s)))ds. 


En supposant que v(s, x) est bien définie aussi pour s € [fo — €, to]. On a donc 


k CS, rase [° [7 v(s', TR) (s— s')ds'ds. 
to 


T+e p(s+e)At 
a J (ST (DO (s-s')dsds’, 
to— 


S—E)Vto 
où a V b = max(a, b) et a A b = min(a, b). Dans la dernière intégrale, nous pouvons 


remplacer s par s’ et en même temps s’ par s. Donc, compte tenu de la relation 


8 s-s)= 861 (s — 5), 


on a 
t+e p(s+e)At fe] ; 
| Se YT (DA (s-s'dsds = 
to—E J(S—-E)V to 
t+e pi(s+e)At fe] fe] ; 
v(S, TS DO ts S)ds'ds. 
to— (S—-E)Vto 
On a donc 


ho Vi; 16 rase f” v(s, DES — s'\ds'ds + Ri(E) + RE), 
lo to+ 


lo+E 
R (€) = Ve [” U(s, ge (sDOT.(s-s')ds'ds, 
to— (s—-E)Vt 


t+e p(s+e)At 
R2(€) = [. Î v(s, VE (SDS s')ds'ds. 
é= S—E 


D'autre part, comme 


S+ 
[ 
jo D ts s'\ds'=1, 
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on a 


[ v(S, y! doas= f° [” v(s, TD ts sds'ds = 
Lo to = 


-f" cf v(s, 1. 1) (s-s)ds'ds + R3(E) + RE), 
to+ SE 


lo+E LA 
R3(€) = L vs, (S)ds, Rate) = Î vs, VIS) ds. 
to l—E 


Des égalités établies ci-dessus, on obtient 


t 
Î Œ(S, ÿ. (s)) — v(s, VE (S))ds = = (3.53) 
lo 


f" cf &(s PES) - vs, TE DE (s-sds'ds+ 
Lo+ S—E 
+Ri(e) + R2(e) — R3(e) + RA(e). 


On va examiner la différence intégrale du second membre de (3.53). On remarque 
d’abord que 


os, rs) - vs, DIS Lys) -rÉ OI. 


D'autre part, on a 
FD a1=| [vf ras) = 


<|s—s'|suplv|<esup|v|. 


Par conséquent, on a 
Lots, y (D — v(s, 71 <ELsuplvl, 


ce qui nous donne 


[ [7 WTA) - LE, TE DE, (s- s)ds'ds| < EL(- ty) suplvl. 
to+E J S—-E 


Imane BAZINE 61 Université 8 Mai 1945-Guelma 


Chapitre 3. Affaiblissement de la condition pour l'estimation dans W1 … 


ou encore 


t 
| Î WEST) vs, 7 (SD ds| < EL(+— to) sup|v|+ R (€) + Ro(e) + Rae) + Ra(®). 
to 


[el 


Donc, la fonction Z(9 = f} 1y1i) 


(s) — y(s)1ds vérifie l'inégalité 
Z'(n) < LZ( + eL(t-— to) sup|v| + Ri(e) + Ri(e) + R3(e) + Rae). 


d'où 


t 
Z(t) < eLsuplel | (s— te CS 4s+ R'(e) 
lo 


où R'(e) = (Ri(e)+ R(e)+ R3(e) + Ra(€)) [ " elU-S) gs ce qui nous permet d’obtenir (3.52). 


Démonstration. du lemme 3.5.2. On a 
17 os Xoï 9 — Yo, oi O1 < 17 (Ko, xo; 0) — {5 Go, Xo: DI+ 17 Go xo: D) — Yo; Xos D 


On commence par l'estimation du premier terme du second membre. Comme v ap- 
partient à L®(0, 1, WL1(Q;R")), il existe une constante L telle que |v(#,x) — v(t,x')| < 


L|x- x!| pour tout (f, x, x’) € [0, r] x Q x Q et ceci est valable aussi pour vlel et donc on a 
fe] re El fee el 
1 Go, 0; D — 7 É (o, 0; D = | Î (sy) LEGS dsl 
to 


t t 
PI ROC OT IC) ES Î ls, 7109) - vs, ps) ds 
to lo 
t t 
<L Î RO PIOIT EX Î @s SD - vs SD ds 
to to 


Par la définition de vl£l on a 


ul (s, (5) - va, 7 (91 < 2EL 
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et donc la fonction Y(#) = f : ly£(s) — *. (s)|ds vérifie l'inégalité 


Y'( < LY(#) +2eEL(t — to), 


t 
Y(n <2eL | (s-toe Vds, 
lo 


ce qui nous permet d'obtenir 


t 
(7 Géo Xoï D — Yi os os I 2EL(t- +L | (ste ds). (8.54 
lo 


Des inégalités (3.52) et (3.54) on parvient à (3.51). 


Lemme 3.5.4. On a 


mes(A,;7) 0 pour Ô —0. (9255) 


Démonstration. Soit (f, x) € A7. Alors par la définition de À; =, il existe un € €]0, Êl et 


[e] 


un point (t*,x*) € V5(Ao) tels que (f*,x*)ey{, 


D’après le lemme 3.5.2 on a 
CE, x) — Yet, (OI < CE 
où C est une constante positive. On considère la caractéristique y(1,», alors on a 
IC, 27) ya HE Ce. 


On pose 


Yan) =", 
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Comme (t*,x*)e Vs(Ao),ona 


(1°, € V5+ce(Ao) 


on a donc 
(£, x) € E(Vi+ce(Ao)) € E(V5: cs Ao)) 
ou 
A5,7 © EUVÿi ce (A0). (3.56) 
La propriété (3.55) résulte de la relation (3.56) et de la condition A. [ 


Proposition 3.5.1. On suppose que 


ve L® (0,00; WL (Q;R")), (3.57) 


loc 


Soit u la solution de l'équation (3.8) avec les conditions (3.9)-(3.10). Alors sous les 


mêmes hypothèses complémentaires du théorème 3.2.1, on a 


ul rw, con < Max(A, sup ess B(s)), (3.58) 
O<s<t 
où 
La 
A = a q Ua oem, nt IV vla om ant Dvl1e (0,10) ” 
al 
q 


x [|| Uollwz, co) +1 MAPTATAONE 


En “ai 
B(s) = Un TUSlzaszwg op IV law opt Dvlio(s Le), 


= GE 
X [241 (S) lu crc) HSE IF zacszwez, ol: 


Démonstration. On remarque que l'inégalité (3.48) ne dépend pas de 6, ce qui nous 


permet de passer à la limite pour Ô — 0, de sorte que, à l’aide du lemme 3.5.4, on ob- 


tiendra (3.58). 
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3.6 Démonstration du théorème 3.2.1 


Rappelons les conditions (3.2)-(3.3), selon lesquelles, quel que soit f>0,ona 
veLI(0,Wi(Q),  V:veL1(0,fW}1(Q)). (3.59) 
On introduit les champs de vecteurs approchés vl"l par 


V(E, x) si [v(£,)luw oo < 7 
DAS (3.60) 


U(E,X) : ; 
ane À Emo 71 


Soit 


(os %0) € (0x DUT U 1xr-H})=z.. 
te[0,r] 


Désignons par y(o, Xo; t) et Yn(fo, Xo; f) la solution de l'équation intégrale 


t 
y = “+ [ v(s,7(s))ds, 


lo 


et 


t 
Yn(®) = Xo +/ vs, yn(s)ds 


lo 
respectivement. 


D’après la proposition 3.5.1 l'équation 


ô 
_ +V-(uvl")=gu+f dans Q (3.61) 


avec les conditions (3.9)-(3.11) admet une solution et une seule, que nous notons ul". 
En outre, le problème (3.8) avec (3.9)-(3.11) lui aussi admet une solution, notée u, qui 


sera construite par la famille d'équations sur les courbes y 


d 
a ETC = —(V:v(,7(0)u(, y) + 8, 7)u(, y) + F7), (3.62) 
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Uo(to, Xo) Si to = 0 
U(to, Xo) = ‘ (3.63) 


u1(fo, Xo) Si t >0 


0) = Yo, Xo; t). 


On a le lemme suivant. 


Lemme 3.6.1. Les solutions ul"! du problème (3.61) avec (3.9)-3.11) convergent 


presque partout dans [0, t] x Q vers la solution u construite par la famille d'équations 


(3.62). 
Pour démontrer le lemme 3.6.1, nous avons besoin du lemme suivant 


Lemme 3.6.2. Pour presque tout (to, Xo) € Z-, les fonctions Yn(to, Xo; t) convergent vers 


la fonction y(to, Xo; t) pour n — co. 
Démonstration. On pose 


Sn) = Yn(to; Xo; ?) — Yo, Xo; t). 


t 
80= | GI (s, Ya(s)) — v(s,y(s)ds. 
lo 


Or,ona 
D, ya) = 0, YONNE (6,20) = 2, O1 + ol, 7) = 0, y). 


On rappelle que, comme v(f, x) est lipschitzienne par rapport à x, v!"l(f, x) sont elles 


aussi lipschitziennes et il existe une fonction L(f) indépendante de n telle que 


Lol (s, x) ll (6, x) < L(bIx-x|  Vx,x'eQ. 
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On a donc 


ul, n(0) - VO TON < LD gn (D. 


D'autre part, de la définition de vl"\(#, x), il résulte immédiatement que 
ul, 70) - ve, 7) < OI max, [2 (4, yo — 
On a donc 


t 
Bento | [L(s)1gn(s)1+17(5)1 max(O, | v(s, Il ça) — Dlds, 
lo 


d’où, de manière usuelle, il résulte que |[g,(#)| tend vers 0. 


Démonstration. du lemme 3.6.1 
Soit ul"] la solution de notre équation. On va montrer que u!"| convergent vers u 


au moins presque partout. L'idée fondamentale est d'estimer 
jull(e, x) ur, x. 


Or, (f, x) est traversé par une caractéristique, disons y = (fo, Xo; t). Donc, on peut écrire 


Qull(r, x) — u(t, x)| dans la forme 


[ul Cr, yo, xo: 0) — u(t,Y(to, Xo: D). 


Qu (6, yo, Xo; 0) — u(t, 7 (to, Xo: 1 < 
< [ul (,7 0, Xo: 9) — ul (£, Yalto, %oi )1 + [ul (6, Ya to, Xo: 0) — u(t, (0, Xo: D) 


On remarque que le terme [ul CE, Yo, Xo; t)) — u(t, (Ko, Xo; t))| peut être majoré 
par 


t 
Î Qc, Dalle, 7,6) ut, )1+ 
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+c(E, Ya) cr, ya, YEDI+ IE, Ya) - fE,yEDDadr, 


en appliquant le lemme de Gronwall, on obtient 
in] , ; filet ynt")idt" 
[UE Yo, Xo; 0) — U(E,7(Lo, Xo; I < e”s d x 


t 
« | Ac ,Yn))- ct, YCDlUE,YEDI+If EYE) FU, y) Dar! 


par conséquent 
ul, x) —u(r, 21 < ul, Ia co lY D a+ 


t ! ! l 
selle lieerdt il (le, cola, TE DIE IS, yen) 7 Idr. 
S 


D’après la proposition 3.5.1, on constate que ||ul”1|, oÆwz,(0) St uniformement bor- 
née. D'autre part, le lemme 3.6.1 garantit la convergence de y, (f) à y(f lorsque ñn — co. 
Donc, on déduit que les solutions ul" convergent presque partout dans [0, f] x Q vers 
la solution u construite par la famille d'équations (3.62). 


De plus, ul" vérifie, pour toute p € €!([0, f] x Q) telle que 


pU,x)=0 sur ( Ü (xr:())u(A x), 


O<t<t 
l'égalité 
n 
im 0P | in] (n] 
dt | [u re À -Vo+(V:v  )p-cp+uop)ldx = 
0  Q 
É 


- | Î u1p(f,x)dtdas. 
0 T_(9 


En passant à la limite pour n — æ,ona 
n n 


far [ue +o.vo+go+ roiax= | Î W-mupdsat- | up(0, max. (3.64) 
0 Q 0 T_() Q 
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Donc, u vérifie l'équation (3.8) et les conditions (3.9)-(3.11) presque partout dans les 


domaines relatifs. 


Démonstration du théorème 3.2.1 D’après la proposition 3.5.1, ul"l vérifie l'estimation 
(3.58). En vertu du lemme 3.6.1 et en passante à la limite lorsque n — co, on trouve 


(3.21). 
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C2 


MANS la présente thèse, nous avons considéré l'équation de transport linéaire 


avec la condition d'entrée. Le résultat présenté dans le chapitre 1 constitue 
une contributions valide à la recherche sur les équations de transport, en particulier 
pour l'élaboration de l'estimation dans L® du gradient de la solution qui joue un rôle 
crucial pour l’étude des équations non linéaires. Quant à la seconde partie, qui est ex- 
plorée dans le deuxième chapitre de cette thèse, est dédiée à l'application du premier 
résultat pour démontrer l'existence et l’unicité de la solution d’un système issu d’un 
modèle du mouvement de l’air avec la transition de phase de l’eau (modèle proposé 
dans [SEL 11]). En outre, dans le troisième chapitre, nous avons amélioré le premier 
résultat présenté dans le premier chapitre. Plus précisémment, nous avons affaibli les 
conditions posées dans le premier résultat pour obtenir un résultat analogue sur l’es- 
timation dans W} de la solution de l’équation de transport. 

Les résultats principaux de ce travail ne sont pas nécessairement optimaux dans le 
sens strict. Mais les nouvelles techniques développées dans notre étude nous semble 
suffisemment efficaces. Nous pouvons donc espérer qu’elles donneront des contribu- 
tions utiles pour les futurs recherches sur les équations de transport avec l'entrée. 

De point de vue d'application, nous avons réussi à résoudre le système d'équations 
non linéaires des densités de l’air sec, de la vapeur d’eau, de l’eau liquide contenue 


dans les gouttelettes, système d'équations avec la condition d'entrée, ce qui constitue 


70 


Conclusion et perspectives 


la généralisation du résultat de la second partie du travail [ASC 14]. Nous pouvons donc 
croire que notre résultat ouvre une voie vers la résolution du système d'équations d’un 
modèle général du mouvement de l’air avec la transition de phase de l’eau proposé 


dans [SEL 11]. 
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